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entrance beam of light and to provide a first reflected beam of 
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intercept the first reflected beam of light from the first deform- 
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and a signal processing and control system configured to 
communicate with the first and second deformable mirrors. 
The first deformable mirror, the second deformable mirror 
and the signal processing and control system together provide 
a large amplitude light modulation range to provide an 
actively shaped optical beam. 
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ACTIVE BEAM SHAPING SYSTEM AND 
METHOD USING SEQUENTIAL 
DEFORMABLE MIRRORS 

FEDERAL FUNDING BY THE U.S. 5 

GOVERNMENT 

This invention was made with Government support of 
Grant No. NAS7-03001, awarded by the National Aeronau- 
tics and Space Administration (NASA). The U.S. Govern- to 
ment has certain rights in this invention. 

BACKGROUND 

1. Field of Invention 15 

The field of the currently claimed embodiments of this 

invention relates to active optical beam shaping systems and 
methods, and more particularly to active optical beam shap- 
ing systems and methods using sequential deformable mir- 
rors. 20 

2. Discussion of Related Art 

Exo-planetary systems that are directly imaged using exist- 
ing facilities (Marois et al. 2008; Kalas et al. 2008; Lagrange 
et al. 2010) give a unique laboratory to constrain planetary 
formation at wide separations (Rafikov 2005; Dodson-Rob- 25 
inson et al. 2009; Kratter et al. 2010; Johnson et al. 2010), to 
study the planetary luminosity distribution at critical young 
ages (Spiegel & Burrows 2012; Fortney et al. 2008) and the 
atmospheric properties of low surface gravity objects (Bar- 
man et al. 201 lb, a; Marley et al. 2010; Madhusudhan et al. 30 
2011). Upcoming surveys, conducted with instruments spe- 
cifically designed for high-contrast (Dohlen et al. 2006; Gra- 
ham et al. 2007; Hinkley et al. 201 1 ), will unravel the bulk of 
this population of self-luminous Jovian planets and provide 
an unprecedented understanding of their formation history'. 35 
Such instruments will reach the contrast required to achieve 
their scientific goals by combining Extreme Adaptive Optics 
systems (Ex -AO, Poyneer & V’eran (2005 )), optimized coro- 
nagraphs (Soummer et al. 2011; Guyon 2003; Rouan et al. 
2000) and nanometer class wavefront calibration (Sauvage et 40 
al. 2007; Wallace et al. 2009; Pueyo et al. 2010). In the future, 
high-contrast instruments on Extremely Large Telescopes 
will focus on probing planetary formation in distant star 
forming regions (Macintosh et al. 2006), characterizing both 
the spectra of cooler gas giants (V’erinaudetal. 2010) and the 45 
reflected light of planets in the habitable zone of low mass 
stars. The formidable contrast necessary to investigate the 
presence of biomarkers at the surface of earth analogs (>1 0 10 ) 
cannot be achieved from the ground beneath atmospheric 
turbulence and will require dedicated space-based instru- 50 
ments (Guyon 2005). 

The coronagraphs that will equip upcoming Ex-AO instru- 
ments on 8 meter class telescopes have been designed for 
contrasts of at most ~10 -7 . Secondary support structure (or 
spiders: 4 struts each 1 cm wide, -0.3% of the total pupil 55 
diameter in the case of Gemini South) have a small impact on 
starlight extinction at such levels of contrasts. In this case, 
coronagraphs have thus been optimized on circularly sym- 
metric apertures, which only take into account the central 
obscuration (Soummer et al. 201 1). However, high-contrast 60 
instrumentation on future observatories will not benefit from 
such gentle circumstances. ELTs will have to support a sub- 
stantially heavier secondary than 8 meter class observatories 
do, and over larger lengths : as a consequence, the relative area 
covered by the secondary support will increase by a factor of 65 
10(30 cm wide spiders, occupying -3% of the pupil diameter 
in the case of TMT). This will degrade the contrast of coro- 
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nagraphs only designed for circularly obscured geometries by 
a factor - 1 00, when the actual envisioned contrast for an ELT 
exo-planet imager can be as low as ~10 -8 (Macintosh et al. 
2006). While the trade-offs associated with minimization of 
spider width in the space-based case have yet to be explored, 
secondary support structures will certainly hamper the con- 
trast depth of coronagraphic instruments of such observato- 
ries at levels that are well above the 1 0 1 0 contrast requirement. 
As a consequence, telescope architectures currently envi- 
sioned for direct characterization of exo-earths consist of 
monolithic, off-axis, and thus im-obscured, telescopes 
(Guyon et al. 2008; Trauger et al. 2010). Coronagraphs for 
such architectures take advantage of the pupil symmetry to 
reach a theoretical contrast of ten orders of magnitude 
(Guyon et al. 2005; Vanderbei et al. 2003a, b; Kasdin et al. 
2005; Mawet et al, 2010; Kuchner & Traub 2002; Soummer et 
al. 2003). However, using obscured on-axis and/or segmented 
apertures take full advantage of the limited real estate asso- 
ciated with a given launch vehicle and can allow larger aper- 
tures that increase the scientific return of space-based direct 
imaging survey. Recent solutions can mitigate the presence of 
secondary support structures in on-axis apertures. However 
these concepts present practical limitations: APLCs on arbi- 
trary apertures (Soummer et al. 2009) and Shaped Pupils 
(Carlotti et al. 201 1 ) suffer from throughput loss for very high 
contrast designs, and PIAAMCM (Guyon et al. 2010a) rely 
on a phase mask technology whose chromatic properties have 
not yet been hilly characterized. Moreover, segmentation will 
further complicate the structure of the telescope’s pupil: both 
the amplitude discontinuities created by the segments gaps 
and the phase discontinuities resulting from imperfect phas- 
ing will thus further degrade coronagraphic contrast. Devis- 
ing a practical solution for broadband coronagraphy on asym- 
metric, unfriendly apertures is an outstanding problem in high 
contrast instrumentation. Therefore, there remains a need for 
improved beam shaping systems and methods. 

SUMMARY 

An active optical beam shaping system according to an 
embodiment of the current invention includes a first deform- 
able mirror arranged to at least partially intercept an entrance 
beam of light and to provide a first reflected beam of light, a 
second defonnable mirror arranged to at least partially inter- 
cept the first reflected beam of light from the first deformable 
mirror and to provide a second reflected beam of light, and a 
signal processing and control system configured to commu- 
nicate with the first and second defonnable mirrors. The 
signal processing and control system is configured to provide 
control signals to the first deformable mirror so as to config- 
ure a reflecting surface of the first deformable mirror to sub- 
stantially conform to a calculated surface shape so that the 
first reflected beam of light will have a predetermined ampli- 
tude profile. The signal processing and control system is 
further configured to provide control signals to the second 
deformable mirror to configure a reflecting surface of the 
second defonnable mirror to substantially conform to a cal- 
culated surface shape so that a phase delay profile across the 
second reflected beam of light will be substantially constant 
when a phase delay profile across the entrance beam of light 
is substantially constant. The first deformable mirror, the 
second defonnable mirror and the signal processing and con- 
trol system together provide a large amplitude light modula- 
tion range to provide an actively shaped optical beam. 

An optical transmitter according to an embodiment of the 
current invention includes an optical source, and an optical 
modulator arranged in an optical path of the optical source. 
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The optical modulator includes a first deformable mirror 
arranged to at least partially intercept an entrance beam of 
light from the optical source and to provide a first reflected 
beam of light, a second deformable mirror arranged to at least 
partially intercept the first reflected beam of light from the 5 
first deformable mirror and to provide a second reflected 
beam of light, and a signal processing and control system 
configured to communicate with the first and second deform- 
able mirrors. The signal processing and control system is 
configured to provide control signals to the first deformable to 
mirror and the second deformable mirror to provide at least a 
portion of the second reflected beam of light as an output 
beam having at least one of a selectable amplitude profile or 
selectable phase profile. 

An optical receiver according to an embodiment of the 15 
current invention includes an optical detector, and an optical 
filter arranged in an optical path of the optical detector. The 
optical filter includes a first deformable mirror arranged to at 
least partially intercept an entrance beam of light being 
detected to provide a first reflected beam of light, a second 20 
deformable mirror arranged to at least partially intercept the 
first reflected beam of light from the first deformable mirror 
and to provide a second reflected beam of light, and a signal 
processing and control system configured to communicate 
with the first and second deformable mirrors. The signal 25 
processing and control system is configured to provide con- 
trol signals to the first deformable mirror and the second 
deformable mirror to provide at least a portion of the second 
reflected beam of light as an output beam directed to the 
optical detector having at least one of a selectable amplitude 30 
profile or selectable phase profile. 

An optical communication system according to an embodi- 
ment of the current invention includes an optical transmitter, 
an optical transmission path optically coupled to the optical 
transmitter, and an optical receiver optically coupled to the 35 
optical transmission path. At least one of the optical transmit- 
ter, the optical transmission path, or the optical receiver 
includes an active optical beam shaper. The active optical 
beam shaper includes a first deformable mirror arranged to at 
least partially intercept an entrance beam of light to provide a 40 
first reflected beam of light, a second deformable mirror 
arranged to at least partially intercept the first reflected beam 
of light from the first deformable mirror and to provide a 
second reflected beam of light, and a signal processing and 
control system configured to communicate with the first and 45 
second deformable mirrors. The signal processing and con- 
trol system is configured to provide control signals to the first 
deformable mirror and the second deformable mirror to pro- 
vide at least a portion of the second reflected beam of light as 
an output beam having at least one of a selectable amplitude 50 
profile or selectable phase profile. 

An optical telescope according to an embodiment of the 
current invention includes a light collection and magnifica- 
tion system, an active optical beam shaper arranged in an 
optical path of light collected and magnified by the light 55 
collection and magnification system, and an optical detection 
system arranged to receive a beam of light from the active 
optical beam shaper. The active optical beam shaper includes 
a first deformable mirror arranged to at least partially inter- 
cept an entrance beam of light to provide a first reflected beam 60 
of light, a second deformable mirror arranged to at least 
partially intercept the first reflected beam of light from the 
first deformable mirror and to provide a second reflected 
beam of light, and a signal processing and control system 
configured to co mmu nicate with the first and second deform- 65 
able mirrors. The signal processing and control system is 
configured to provide control signals to the first deformable 
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mirror and the second deformable mirror to provide at least a 
portion of the second reflected beam of light as an output 
beam having at least one of a selectable amplitude profile or 
selectable phase profile. 

A solar furnace according to an embodiment of the current 
invention includes a light collection system, an active optical 
beam shaper arranged in an optical path of light collected by 
the light collection system, and an optical focusing system 
arranged to receive a beam of light from the active optical 
beam shaper to be focused on an object to be heated. The 
active optical beam shaper includes a first deformable mirror 
arranged to at least partially intercept an entrance beam of 
light to provide a first reflected beam of light, a second 
deformable mirror arranged to at least partially intercept the 
first reflected beam of light from the first deformable mirror 
and to provide a second reflected beam of light, and a signal 
processing and control system configured to communicate 
with the first and second deformable mirrors. The signal 
processing and control system is configured to provide con- 
trol signals to the first deformable mirror and the second 
deformable mirror to provide at least a portion of the second 
reflected beam of light as an output beam having at least one 
of a selectable amplitude profile or selectable phase profile. 

An optical pulse-shaping system according to an embodi- 
ment of the current invention includes an optical pulse source, 
a first diffraction grating disposed in an optical path of the 
optical pulse source, an active optical beam shaper disposed 
in an optical path of light diffracted from the first diffraction 
grating, and a second diffraction grating disposed in an opti- 
cal path of light output from the active optical beam shaper. 
The active optical beam shaper includes a first deformable 
mirror arranged to at least partially intercept an entrance 
beam of light to provide a first reflected beam of light, a 
second deformable mirror arranged to at least partially inter- 
cept the first reflected beam of light from the first deformable 
mirror and to provide a second reflected beam of light, and a 
signal processing and control system configured to commu- 
nicate with the first and second deformable mirrors. The 
signal processing and control system is configured to provide 
control signals to the first deformable mirror and the second 
deformable mirror to provide at least a portion of the second 
reflected beam of light as an output beam having at least one 
of a selectable amplitude profile or selectable phase profile. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Further objectives and advantages will become apparent 
from a consideration of the description, drawings, and 
examples. 

FIG. 1A is a schematic illustration of an active optical 
beam shaping system according to an embodiment of the 
current invention. 

FIG. IB is a schematic illustration of an optical transmitter 
according to an embodiment of the current invention. 

FIG. 1C is a schematic illustration of laser beam launch 
system according to an embodiment of the current invention. 

FIG. ID is a schematic illustration of an optical receiver 
according to an embodiment of the current invention. 

FIG. IE is a schematic illustration of a solar furnace 
according to an embodiment of the current invention. 

FIG. IF is a schematic of a beam shaping device according 
to an embodiment of the current invention. Some advantages 
of beam shaping can include: (1) Can put more energy in a 
pulse. These can be ultra-short, high energy femtosecond 
pulses. Can include applications to laser fusion. (2) Can 
actively correct for beam aberrations that do not degrade the 
gratings. 
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FIG. 1G is a schematic illustration of a telescope according 
to an embodiment of the current invention. Top: Architecture 
applicable to future exo-earth imaging missions: a monolithic 
un-obscured telescope feeds a coronagraph designed on a 
circular aperture. The wavefront errors are corrected using 
two sequential Deformable Mirrors (DMs) that are controlled 
using a quasi-linear feedback loop based on image plane 
diagnostics. The propagation between optical surfaces, 
between the DMs in particular is assumed to occur in the 
Fresnel regime. Bottom: ACAD solution: the coronagraph is 
designed for a circular geometry around the central obscura- 
tion. The two sequential DMs are controlled in the non-linear 
regime based on a pupil plane cost function. The propagation 
between the DMs is modeled using ray optics, and we con- 
duct a quantitative one-dimensional analysis of the diffrac- 
tion artifacts. 

FIG. 1H is a schematic illustration to help explain some 
concepts of some embodiments of the current invention. 

FIG. 2 shows optimally designed Apodized Pupil Lyot 
Coronagraphs on circularly obscured apertures. With fixed 
obscuration ratio, size of opaque focal plane mask, IWA and 
OWA, our linear programming approach yields solutions 
with theoretical contrast below 10 10 . All PSFs shown on this 
figure are monochromatic for \='k 0 . When all other quantities 
remain equal and the central obscuration ratio increases (from 
10% in the top panel to 20% in the middle panel), then the 
solution becomes more oscillatory (e.g less feasible) and the 
contrast constraint has to be relaxed. Eventually the optimizer 
does not find a solution and the size of the opaque focal plane 
mask (and thus the IWA) has to be increased (central obscu- 
ration of 30% in the bottom panel). On the right hand side, we 
present our results in two configurations: when the apodiza- 
tion is achieved using at grayscale screen (APLC), “on-sky” 
Z/D bottom x-axis, and when the apodization is achieved via 
two pupil remapping mirrors (PIAAC), “on-sky” Z/D top 
x-axis. We adopt this presentation to show that ACAD is 
“coronagraph independent” and that it can be applied to coro- 
nagraphs with high throughput and small IWA. 

FIG. 3 shows a schematic of the optical system considered: 
the telescope aperture is followed by two sequential Deform- 
able Mirrors (DMs) in non-conjugate planes whose purpose 
is to remap the pupil discontinuities. The beam then enters a 
coronagraph to suppress the bulk of the starlight: in this figure 
we show an Apodized Lyot Pupil Coronagraph (APLC). This 
is the coronagraphic architecture we consider for the remain- 
der of the paper but we stress that the method presented herein 
is applicable to any coronagraph. 

FIG. 4 Top: In the ideal case, the two DMs would hilly 
remap all the discontinuities in the telescope’s aperture to 
feed a hilly uniform beam to the coronagraph. However this 
would require discontinuities in the mirror’s curvatures 
which are not achieved in practice. Moreover solving the 
Monge-Ampere Equation in this direction is a difficult exer- 
cise as the right hand side of Eq. 26 presents an implicit 
dependence on the solution h 2 . 

Bottom: We circumvent this problem by solving the 
reverse problem, where the the input beam is now uniform 
and the implicit dependence drops out. Moreover we taper the 
edges of the discontinuities by convolving the target field 
A(x 15 y ,) by a gaussian of full width at half maximum w 
(w=50 cycles per aperture in this figure). 

FIG. 5 shows forward and reverse coordinate systems, with 
their respective forward and reverse coordinate transforms, in 
the case of the folded system studied in this example. Note 
that this figure shows normalized units. As explained in the 
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body of the text the correspondence between normalized and 
real units scales as follow: (x,,y l )=(DX I ,DY I ), (f„g,)=(DF„ 
DG,), 

5 



10 where i=l, 2, D is the aperture diameter and Z the separation 
between DMs. We solve the Monge-Ampere Equation in 
normalized coordinates and then apply the scalings in order to 
find the true DM shapes. 

FIG. 6 shows the virtual field at Ml when solving the 
15 reverse problem, A(X 1 ,Y 1 ) is the desired apodization (top), 
A„(X 1 ,Y 1 ) is the apodization obtained after solving the 
Monge-Ampere Equation (center). The bottom panel shows 
the difference between the two quantities: the bulk of the 
energy in the residual is located in high spatial frequencies 
20 that cannot be controlled by the DMs. 

FIG. 7 shows boundary conditions seen in the horizontal 
remapped space. Should the boundary conditions have 
strictly been enforced by our solver then 

30 

The remapping function obtained with our solutions is very 
close to these theoretical boundary conditions and the residu- 
als can easily be mitigated by sacrificing the edge rows and 
columns of actuators on each DM. 

FIG. 8 shows boundary conditions seen in the vertical 
remapped space. Should the boundary conditions have 
strictly been enforced by our solver then 



45 

The remapping function obtained with our solutions is very 
close to these theoretical boundary conditions and the residu- 
als can easily be mitigated by sacrificing the edge rows and 
columns of actuators on each DM. 

50 FIG. 9 shows beam amplitude before and after the DMs in 
the case of a geometry similar to JWST. We chose to solve the 
reverse problem over a square using a Fourier basis set and not 
strictly enforcing boundary conditions. This results in small 
distortions of the edges of the actual aperture in the vicinity of 
55 segments gaps and secondary supports. We address this prob- 
lem by slightly oversizing the secondary obscuration and 
undersizing the aperture edge in the coronagraph. 

FIG. 10, Left: Comparison of edge diffraction between 
PIAA and Fresnel propagation. Because the discontinuities in 
60 the pupil occur at the location where T^cl , a PIAA amplifies 
the chromatic ringing when compared to a more classical 
propagation. Right: Comparison of edge diffraction between 
ACAD and Fresnel propagation. Because the discontinuities 
in the pupil occur at the location where T X >1 , a ACAD damps 
65 the chromatic ringing when compared to a more classical 
propagation. These simulations were carried out for a one 
dimensional aperture. In this case the Monge-Ampere Equa- 
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tion can be solved using finite elements and the calculation of 
the diffractive effects reduces to the evaluation of various 
Fresnel special functions at varying wavelength. The full 
two-dimensional problem is not separable and requires the 
development of novel numerical tools. 5 

FIG. 11 Results obtained when applying our approach to a 
geometry similar to JWST. We used two 3 cm DMs of 64 
actuators separated by 1 m. Their maximal surface deforma- 
tion is 1.1 pm, well within the stroke limit of current DM 
technologies. The residual light in the corrected PSF follows to 
the secondary support structures and can potentially be fur- 
ther cancelled by controlling the DMs using an image plane 
based cost function, see FIG. 23. 

FIG. 12 Case of JWST: Radial average obtained when 
applying ACAD. We used two 3 cm DMs of 64 actuators 15 
separated by 1 m. Their maximal surface deformation is 1 . 1 
pm, well within the stroke limit of current DM technologies. 
ACAD yields a gain in contrast of two orders of magnitude, 
and provides contrasts levels similar to upcoming Ex -AO 
instruments, which are designed on much friendlier apertures 20 
geometries. Since ACAD removes the bulk of the light dif- 
fracted by the asymmetric aperture discontinuities, the final 
contrast can be improved by controlling the DMs using and 
image plane based cost function, see FIG. 23. 

FIG. 13 shows results obtained when applying our 25 
approach to a TMT geometry. We used two 3 cm DMs of 64 
actuators separated by 1 m. Their maximal surface deforma- 
tion is 0.9p 0.9 pm, well within the stroke limit of current DM 
technologies. The final contrast is below 10 7 10 7 , in a regime 
favorable for direct imaging of exo-planets with ELTs. 30 

FIG. 14 Case of TMT: radial average obtained when apply- 
ing ACAD. We used two 3 cm DMs of 64 actuators separated 
by 1 m. Their maximal surface deformation is 0.9 p 0.9 pm, 
well within the stroke limit of current DM technologies. The 
final contrast is below 1.0 7 10 7 , in a regime favorable for 35 
direct imaging of exo-planets with ELTs. Since ACAD 
removes the bulk of the light diffracted by the asymmetric 
aperture discontinuities, the final contrast can be enhanced by 
controlling the DMs using and image plane based metric. 

FIG. 15 shows PSFs resulting from ACAD when varying 40 
the number and thickness of secondary support structures 
while maintaining their covered surface constant. The surface 
area covered in this example is 50% greater than in the TMT 
example shown on FIG. 13. As the spiders get thinner their 
impact on raw contrast becomes smaller and the starlight 45 
suppression after DM correction becomes bigger. For a rela- 
tively small number of spiders (<12) the contrast improve- 
ment on each single structure is the dominant phenomenon, 
regardless of the number of spiders. ELTs designed with a 
moderate to large number of thin secondary support structure 50 
(6 to 12) present aperture discontinuities which are easy to 
correct with ACAD. 

FIG. 16 shows radial PSF profiles resulting front ACAD 
when varying the number and thickness of secondary support 
structures while maintaining their covered surface constant. 55 
The surface area covered in this example is 50% greater than 
in the TMT example shown on FIG. 13. As the spiders get 
thinner their impact on raw contrast becomes lesser and the 
starlight suppression after DM correction becomes greater. In 
the 12 spiders example, at latge separations, the average 60 
contrast is an order of magnitude higher than reported on FIG. 

14 . 

FIG. 17 shows PSFs resulting from ACAD when varying 
the number and thickness of secondary support structures 
while maintaining their covered surface constant. The surface 65 
area covered in this example is 50% greater than in the TMT 
example shown on FIG. 13. When the number of spiders 
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increases, they produce a sharp circular diffraction feature at 
N SpjJer Jn X/D. If this number is greater than the size of the 
focal plane mask this structure appears in the high contrast 
zone and is very difficult to correct with ACAD. The bright- 
ness of this structure is mitigated by the fact that the spiders 
are very thin. 

FIG. 18 shows radial PSF profiles resulting from ACAD 
when varying the number and thickness of secondary support 
structures while maintaining their covered surface constant. 
The surface area covered in this example is 50% greater than 
in the TMT example shown on FIG. 13. With a large number 
of spiders the bright ring in the PSF structure located at 
^spiderJ 71 is difficult to correct with ACAD. However, 
since the spiders become thinner their net effect on contrast 
after ACAD remains small. 

FIG. 19 shows PSFs resulting from ACAD when varying 
the number and thickness of secondary support structures. As 
the spiders get thinner their impact on raw contrast becomes 
lesser and the starlight suppression after DM correction 
becomes greater. In this case co was optimized on a very fine 
grid and the aperture we clocked in a favorable direction with 
respect to the Fourier basis. 

FIG. 20 shows PFSs resulting from ACAD when varying 
the number and thickness of secondary support structures. As 
the spiders get thinner their impact on raw contrast becomes 
lesser and the starlight suppression after DM correction 
becomes greater. In this case co was optimized on a very fine 
grid and the aperture we clocked in a favorable direction with 
respect to the Fourier basis. Even for spiders as thick as 0.5% 
of the telescope aperture the designed contrast of the corona- 
graph is retrieved. 

FIG. 21 shows off-axis PSF after ACAD in the case of a 
geometry similar to JWST. The aspheric surface of the DMs 
introduces a slight field-dependent distortion. However the 
core of the PSF is still concentrated within the central airy 
disk and the DMs only have an effect on the PSF tail. Field 
distortion does not thus hamper the detectability of faint 
off-axis sources. 

FIG. 22 shows broadband wavefront correction (20% 
bandwidth around 700 mn) with a single DM in segmented 
telescope with discontinuous surface errors. Top Left: wave- 
front before the coronagraph. Top Right: broadband aberrated 
PSF with DM at rest. Bottom Left: DM surface resulting from 
the wavefront control algorithm. Bottom Right: broadband 
corrected PSF. Note that the wavefront control algorithm 
seeks to compensate for the diffractive artifacts associated 
with the secondary support structures: it attenuates them on 
the right side of the PSF while it strengthens them on the left 
side of the PSF. As a result the DM surface becomes too large 
at the pupil spider’s location and the quasi-linear wavefront 
control algorithm eventually diverges. 

FIG. 23 shows broadband wavefront correction (20% 
bandwidth around 700 mn) in a segmented telescope whose 
pupil has been re-arranged using ACAD. The surface of the 
first DM is set according to the ACAD equations. The surface 
of the second DM is the sum of the ACAD solution and a 
small perturbation calculated using a quasi-linear wavefront 
control algorithm. Top Left: wavefront before the corona- 
graph. Note that the ACAD remapping has compressed the 
wavefront errors near the struts and the segment gaps. Top 
Right: broadband aberrated PSF with DMs set to the ACAD 
solution. Bottom Left: perturbation of DM2’s surface result- 
ing from the wavefront control algorithm. Bottom Right: 
broadband corrected PSF. The wavefront control algorithm 
now yields a DM surface that does not feature prominent 
deformations at the location of the spiders. Most of the DM 
stroke is located at the edge of the segments, at location of the 
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wavefront discontinuities. There, the DM surface eventually 
becomes too large and the quasi-linear wavefront control 
algorithm diverges. However this results in higher contrasts 
than in the absence of ACAD. 

FIG. 24 shows radial average in the half dark plane of the 
PSFs on FIG. 22 and FIG. 23. In the presence of wavefront 
discontinuities corrected using a continuous membrane DM, 
ACAD still yields, over a 20% bandwidth around 700 nm, 
PSF with a contrast 100 larger than in a classical segmented 
telescope. Moreover this figure illustrates that since it is based 
on a true image plane metric, the wavefront control algorithm 
can be used (within the limits of its linear regime) to improve 
upon the ACAD DM shapes derived solving the Monge- 
Ampere Equation. 

FIG. 25 shows further embodiments for higher contrasts 
with ACAD. The blue and orange colors respectively repre- 
sent the current state of the art in wavefront control and the 
work described in the present manuscript, as in FIG. 1G. In 
brown are listed avenues to further the contrasts presented 
herein: 1) combining ACAD with coronagraphs designed on 
segmented and/or on-axis apertures, 2) using diffractive mod- 
els to close a quasi-linear focal plane based loop using a 
metric whose starting point corresponds to the DM shapes 
calculated in the non-linear regime. 

FIG. 26 shows simulation of ACAD in the SR-Fresnel 
approximation. We simulated the propagation between DM1 
and DM2 according to Eq. (59), and simulated the rest of the 
coronagraphic train using the Fresnel approximation. Fresnel 
propagation simply replicates the input pupil’s (top panels) 
deep and sharp discontinuities in the final reimaged pupil 
(right bottom panel). ACAD SR-Fresnel instead transforms 
them into shallow and smooth mid-spatial frequency ripples 
(left bottom panel). The chromatic nature of these oscillations 
sets the ultimate bandwidth of ACAD, but unlike in PIAA 
coronagraphs the beam reshaping here does not amplify these 
oscillations any beyond normal edge ringing for out of pupil 
optics in the Fresnel regime (discussed in Pueyo and Kasdin, 
2007). The plotted slices across the secondary support struc- 
tures (middle panels) show that ripples in the SR-Fresnel 
approximation are neither bigger nor more chromatic than 
hard-edge induced ripples in the Fresnel approximation. This 
enables a broad spectral bandwidth. Note that in order to 
better illustrate the oscillations’ chromatic nature this figure 
has been generated with a large separation between the DMs, 
Z=2 m. For more realistic separations ACAD exhibits a 
weaker chromatic behavior, as illustrated by FIG. 2. Here we 
intentionally omit the coronagraphic mask and Lyot stop to 
highlight that ACAD can operate with any coronagraph. 

FIG. 27 shows ACAD in the SR-Fresnel diffractive optics 
regime +RA-VYC}: We demonstrate a simulated 30% band- 
width contrast of ~3x 1 0 9 and an inner working angle of 2 A/D 
on the AFTA pupil using ACAD with two 1 inch DMs sepa- 
rated by 0.5 m. It demonstrates the weak chromaticity of 
ACAD and its ability to operate with a high throughput, small 
IWA coronagraph on an obscured aperture. 

DETAILED DESCRIPTION 

Some embodiments of the current invention are discussed 
in detail below. In describing embodiments, specific termi- 
nology is employed for the sake of clarity. However, the 
invention is not intended to be limited to the specific termi- 
nology so selected. A person skilled in the relevant art will 
recognize that other equivalent components can be employed 
and other methods developed without departing from the 
broad concepts of the current invention. All references cited 
anywhere in this specification, including the Background and 
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Detailed Description sections, are incorporated by reference 
as if each had been individually incorporated. 

The terms “light” and “optical” are intended to have broad 
meanings to include both visible and non-visible regions of 
5 the electromagnetic spectrum. Visible, infrared, ultraviolet 
and x-ray regions of the electromagnetic spectrum are all 
intended to be included within the definition of the term 
“light”. 

FIG. 1A is a schematic illustration of an active optical 
to beam shaping system 1 00 according to an embodiment of the 
current invention. The optical beam shaping system 100 
includes a first deformable mirror 102 arranged to at least 
partially intercept an entrance beam of light 104 and to pro- 
vide a first reflected beam of light 106, a second deformable 
15 mirror 108 arranged to at least partially intercept the first 
reflected beam of light 106 from the first deformable mirror 
102 and to provide a second reflected beam of light 110, and 
a signal processing and control system 112 configured to 
communicate with the first and second deformable mirrors 
20 (102, 108). The signal processing and control system 112 is 
configured to provide control signals to the first deformable 
mirror 102 so as to configure a reflecting surface of the first 
deformable mirror 102 to substantially conform to a calcu- 
lated surface shape so that the first reflected beam of light 106 
25 will have a predetermined amplitude profile. The signal pro- 
cessing and control system 112 is further configured to pro- 
vide control signals to the second deformable mirror 108 to 
configure a reflecting surface of the second deformable mir- 
ror 108 to substantially conform to a calculated surface shape 
30 so that a phase delay profile across the second reflected beam 
of light 110 will be substantially constant when a phase delay 
profile across the entrance beam of light 104 is substantially 
constant. When the phase delay of entrance beam of light 104 
is not constant the signal processing and control system 112 
35 can control the second deformable mirror 108 to make the 
phase delay of the second reflected beam of light 110 con- 
stant. The first deformable mirror 1 02, the second deformable 
mirror 108 and the signal processing and control system 112 
together provide a large amplitude light modulation range to 
40 provide an actively shaped optical beam. 

The term “deformable mirror” refers to an assembly that 
includes a reflecting surface as well as structural components 
that effect reconfiguration of the surface to substantially 
match calculated surface shapes. For example, the reflecting 
45 surface can be reconfigured into a plurality of complex non- 
spherical, non-parabolic and non-ellipsoidal surface shapes 
that can include a plurality of local convex and/or concave 
localized regions. The structural components of the deform- 
able mirrors can include, but are not limited to, actuators. 

50 The signal processing and control system 112 can be an 
application-specific and/or a general programmable system. 
For example, it can be implemented on a computer system 
and/or implemented on dedicated hardware such as an ASIC 
and/or FPGA. The signal processing system and control 112 
55 can also include memory and/or data storage in addition to 
data processing components, such as, but not limited to a CPU 
and/or GPU. The signal processing and control functions can 
be incorporated together or performed on separate compo- 
nents, for example. The communication paths with the 
60 deformable mirrors can be hard wired, such as electrically 
wired or fiber optic, and/or wireless such as radio frequency 
(RF) and/or optical communications links, for example. 

In an embodiment of the current invention, the large ampli- 
tude modulation range can be at least 1% of a fully illumi- 
65 nated portion of the entrance beam 1 04 of light. In an embodi- 
ment of the current invention, the large amplitude modulation 
range can be substantially an entire range to permit selectable 
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amplitude modulation of portions of the entrance beam 1 04 of 
light anywhere from substantially fully attenuated to substan- 
tially hilly illuminated relative to fully illuminated portions of 
the entrance beam of light 104 . If the spatial scale amplitude 
or phase structures of interest of the incoming beam is (8x,8y ) 5 
-O (8) (where 8 is that distance between actuators, OR size of 
discontinuities between telescope pupil, OR solid angle of a 
fiber when coupling to a fiber, or scale of beam/over what info 
is encoded for laser can, or scale of beam perturbations), then 
the large amplitude modulation regime over this characteris- 10 
tic scale can be defined as follows: 

Define 


12 

Assuming that the characteristic scale of the input beam is 
D 

aU 

(actuator pitch) and that n is “as large as possible” in order not 
to lose generality, eg, 

Nad 


dh 

dr 


the focal gradient of both mirror’s surface at any arbitrary 
point on either the input or the output beam (FIG. 1H). 

Define 0^ as the solid angle subtended by the disk or radius 
8 and the distance z between the DMs. 

The radius 8 is a characteristic scale of the beam that is 
sought to be reshaped, Z is a distance between the deformable 
mirrors, defined by a) the two points for which the optical 
surfaces of the DMs are the farthest one to another in the 
horizontal direction, and b) the length between the two tan- 
gent planes to the DM surfaces at those points, lq and h 2 are 
the surface deflections of each DM measured as the relative 
distance from each respective DM surface at the tangent 
planes at maximum horizontal length defined above. 

Then the long amplitude modulation regime is defined 
when significant beam remapping occurs. 


dh 

dr 


* <w 


Eqn (1) 


This rigorous expression can be translated in terms of the 
largest amplitude modulation (“a”) in a simple case. 

Assume we are seeking a ID amplitude modulation: 


1 + flCos 



Then under the assumption a«l one can seek: 


h\{x) = —hi{x) - eCos 



with 


n 2 D 2 
e ~ 4 7T 2 2 ° 

When in the “large amplitude modulation” regime the 
“a«l” breaks down and, this simple solution is not valid, we 
have to solve the Monge-Ampere equation as per the present 
description. 


15 Then folding this into Eqn (1) yields: 



20 

Namely, when N=32 we find a crir =l .2% as previously 
stated in this specification, the boundary between “small” and 
“large” amplitude modulation is around 1%, but the above 
also provides a rigorous derivation. This is the rigorous defi- 
nition that defines the boundary between small and large 
amplitude modulation regime for ALL applications. 

In an embodiment of the current invention, the active opti- 
cal beam shaping system 100 can further include a beam 
30 analyzer 114 arranged to receive at least a portion of the exit 
beam 110 of light reflected from the second deformable mir- 
ror 108 . In alternative embodiments, the beam analyzer 114 
can receive external information regarding the medium or 
environment in which the second beam of light 110 will 
35 propagate. The beam analyzer 114 can be configured to com- 
municate with the signal processing and control system 112 
to provide feedback control of at least one of the first and 
second deformable mirrors ( 102 , 108 ). The beam analyzer 
114 can be hard wired by electrical and/or fiber optic connec- 
40 tions to the signal processing and control system 112 , and/or 
in wireless communication such as by RF and/or optical 
wireless communications links. In some embodiments, the 
beam analyzer 114 can be further configured to estimate: (1) 
the amplitude profile of the beam and/or (2) the phase of the 
45 beam. 

In an embodiment of the current invention, input 116 for a 
desired output beam shape can be provided to the signal 
processing and control system 112 . The input 116 informa- 
tion for the desired output beam shape can be used in con- 
50 junction with feedback signals from said beam analyzer 114 
in some embodiments of the current invention. However, 
embodiments in which no feedback information is used are 
intended to be included within the general concepts of the 
current invention. 

55 In an embodiment of the current invention, the signal pro- 
cessing and control system 112 can be further configured to 
calculate the surface shape of at least one of the first and 
second deformable mirrors ( 102 , 108 ) based on reverse ray 
tracing calculations. 

60 In an embodiment of the current invention, the signal pro- 
cessing and control system 112 can be further configured to 
calculate the surface shape of at least one of the first and 
second deformable mirrors based on a solution to Monge- 
Ampere equations corresponding to the active optical beam 
65 shaping system. 

Further embodiments of the current invention are directed 
to numerous devices and systems that incorporate at least one 
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active optical beam shaping system according to an embodi- 
ment of the current invention. The following will describe 
some such devices and systems; however, one should recog- 
nize that further devices and systems, and modifications 
thereof, can be constructed with one or more active optical 
beam shaping systems according to an embodiment of the 
current invention and are considered to be included within the 
general scope of the current invention. 

FIG. IB is a schematic illustration of an optical transmitter 
200 according to an embodiment of the current invention. The 
optical transmitter 200 includes an optical source 202 and an 
optical modulator 204 arranged in an optical path 206 of the 
optical source 202. The optical modulator 204 can be similar 
to or the same as the active optical beam shaping system 100, 
for example. The optical modulator 204 includes, similar to 
FIG. 1A and not shown in FIG. IB, a first deformable mirror 
arranged to at least partially intercept an entrance beam of 
light from the optical source 202 and to provide a first 
reflected beam of light, a second deformable mirror arranged 
to at least partially intercept the first reflected beam of light 
from said first deformable mirror and to provide a second 
reflected beam of light, and a signal processing and control 
system configured to communicate with the first and second 
deformable mirrors. The signal processing and control sys- 
tem is configured to provide control signals to the first 
deformable mirror and the second deformable mirror to pro- 
vide at least a portion of the second reflected beam of light as 
an output beam 207 having at least one of a selectable ampli- 
hide profile or selectable phase profile. 

The optical transmitter 200 can further include a beam 
analyzer 208 configured to sample a portion of the output 
beam 207. For example, the optical transmitter 200 can 
include a beam splitter 210, such as a partially silvered mirror 
or a dichroic mirror arranged in the output beam 207 to 
redirect a portion of the output beam 207 to beam analyzer 
208. Flowever, the general concepts of the current invention 
are not limited to this particular example. The beam analyzer 
can be, but is not limited to, a camera and/or a wavefront 
sensor, for example. 

The optical transmitter 200 can be used in, but is not limited 
to, optical communications systems for example. Such opti- 
cal communications systems can be free space systems, or 
can be waveguide systems, for example. Waveguide systems 
can include, but are not limited to, fiber optic systems. The 
term “free space” is intended to have a general meaning to 
include vacuum propagation, such as between satellites and/ 
or spacecraft, or through a medium such as the atmosphere, 
and/or water. 

In some embodiments, the signal processing and control 
system of the optical transmitter 200 can be configured to 
provide control signals to the first deformable mirror and the 
second deformable mirror such that the output beam has both 
a selectable amplitude profile and selectable phase profile. In 
some embodiments, the signal processing and control system 
of the optical transmitter 200 can be configured to provide 
control signals to the first deformable mirror and the second 
deformable mirror such that the output beam of light 207 is 
soliton pulses. The soliton boundary conditions can be input 
to the signal processing and control system, for example. For 
a given fiber we do not calculate the soliton shape here. 
However, we can couple to the fiber any prescribed solitons 
and keep accurately their profile in the presence of distur- 
bances. 

The solitons can be coupled into an optical fiber, for 
example, for long distance fiber optic communication without 
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the need for repeaters for signal regeneration, for example. 
This can be useful for transoceanic fiber optic systems, for 
example. 

FIG. 1C is a schematic illustration of laser beam launch 
5 system 220 according to an embodiment of the current inven- 
tion. It includes an optical transmitter 222 and reflecting 
telescope 224. The optical transmitter 222 can be the same or 
similar to optical transmitter 200, for example. Some appli- 
cations of laser beam launch system 220 can include, but are 
not limited to, tactical and strategic shooting down of satel- 
10 lites. 

FIG. ID is a schematic illustration of an optical receiver 
300, according to an embodiment of the current invention. 
The optical receiver 300 includes an optical detector 302 and 
an optical filter 304 arranged in an optical path of the optical 
15 detector 320. The optical filter can be similar to or the same as 
the active optical beam shaping system 100, for example. The 
optical filter 304 includes a first deformable mirror arranged 
to at least partially intercept an entrance beam of light being 
detected to provide a first reflected beam of light, a second 
20 deformable mirror arranged to at least partially intercept the 
first reflected beam of light from the first deformable mirror 
and to provide a second reflected beam of light, and a signal 
processing and control system configured to communicate 
with the first and second deformable mirrors, (Also, see, e.g., 
25 FIG. 1A.) The signal processing and control system is con- 
figured to provide control signals to the first defonnable mir- 
ror and the second defonnable mirror to provide at least a 
portion of the second reflected beam of light as an output 
beam directed to the optical detector to have at least one of a 
30 selectable amplitude profile or selectable phase profile. 

The optical receiver 300 can further include a beam ana- 
lyzer 306 configured to communicate with the optical filter 
304. The beam analyzer 306 can be, but is not limited to, a 
camera and/or wavefront detector, for example. The optical 
35 receiver 300 can further include a beam splitter 308, such as, 
but not limited to, a partially silvered mirror or a dichroic 
mirror to direct a portion of the received beam of light to the 
beam analyzer. The optical receiver 300 can be used in, but is 
not limited to, an optical communication system. The optical 
40 communication system can be a waveguide system and/or a 
free space system. 

In free space communication systems according to some 
embodiments of the current invention, information can be 
encoded in the shape of the beam without loss of light, such as 
45 would occur with temporal modulation. The broad bandwidth 
can be advantageous for communications. 

FIG. IE is a schematic illustration of a solar furnace 400 
according to an embodiment of the current invention. The 
solar furnace 400 includes a light collection system 402, an 
50 active optical beam shaper 404 arranged in an optical path of 
light collected by the light collection system 402, and an 
optical focusing system 406 arranged to receive a beam of 
light from the active optical beam shaper 404 to be focused on 
an object 408 to be heated. The active optical beam shaper 404 
55 can be the same as, or similar to, active optical beam shaping 
system 100. The efficiency is boosted by shaping the solar 
beam precisely to the heat element with losing energy that can 
scatter and also cause damage. 

The following examples describe some embodiments and 
60 some applications in more detail. However, the broad con- 
cepts of the current invention are not limited to the particular 
examples. 

EXAMPLES 

65 

The examples here describe a family of practical applica- 
tions according to some embodiments of the current inven- 
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tion. As their ultimate performances depend strongly on the 
pupil structure, we limit the scope here to a few characteristic 
examples. Full optimization for specific telescope geom- 
etries, as well as other applications, can be conducted as 
needed. 5 

1 Introduction 

Some embodiments of the current invention can take 
advantage of state-of-the-art Deformable Mirrors (DM) in 
modern high-contrast instruments to address the problem of 
pupil amplitude discontinuities for on-axis and/or segmented to 
telescopes. Indeed, coronagraphs are not sufficient to reach 
the high contrast required to image faint exo-planets: wave- 
front control is needed to remove the light scattered by small 
imperfections on the optical surfaces (Brown & Burrows 
1990). Over the past several years, significant progress has 15 
been made in this area, both in the development of new 
algorithms (Borde & Traub 2006; Give’on et al. 2007) and in 
the experimental demonstration of high-contrast imaging 
with a variety of coronagraphs (Give’on et al. 2007; Trauger 
& Traub 2007; Guyon et al. 2010b; Belikov et al, 2011). 20 
These experiments rely on a system with a single Deformable 
Mirror which is controlled based on diagnostics downstream 
of the coronagraph, either at the science camera or as close as 
possible to the end detector (Wallace et al. 2009; Pueyo et al. 
2010). Such configurations are well suited to correct phase 25 
wavefront errors arising from surface roughness but have 
limitations in the presence of pure amplitude errors (reflec- 
tivity), or phase-induced amplitude errors, which result from 
the propagation of surface errors in optics that are not conju- 
gate to the telescope pupil (Shaklan & Green 2006; Pueyo & 30 
Kasdin 2007). Indeed a single DM can only mimic half of the 
spatial frequency content of amplitude errors and compensate 
for them only on one half of the image plane (thus limiting the 
scientific field of view) over a moderate bandwidth. In theory, 
architectures with two sequential Deformable Mirrors, can 35 
circumvent this problem and create a symmetric broadband 
high contrast PSF (Shaklan & Green 2006; Pueyo & Kasdin 
2007). The first demonstration of symmetric dark hole was 
reported in Pueyo et al. (2009) and has since been generalized 
to broadband by Groff et al. (2011). In such experiments the 40 
coronagraph has been designed over a full circular aperture, 
the DM control strategy is based on a linearization of the 
relationship between surface deformations and electrical field 
at the science camera, and the modeling tools underlying the 
control loop consist of classic Fourier and Fresnel propaga- 45 
tors. This is illustrated on the left panel of FIG. 1. As a 
consequence, a wavefront control system composed of two 
sequential Deformable Mirrors is currently the baseline 
architecture of currently envisioned coronagraphic space- 
based instruments (Shaklan et al. 2006; Krist et al. 2011) and 50 
ELT planet imagers (Macintosh et al. 2006). One can thus 
naturally be motivated to investigate if such wavefront control 
systems can be used to cancel the light diffracted by second- 
ary supports and segments in large telescopes, since such 
structures are amplitude errors, albeit large amplitude errors. 55 

The purpose of the current example is to demonstrate that 
indeed a two Deformable Mirror wavefront control system 
can mitigate the impact of the pupil asymmetries, such as 
spiders and segments, on contrast and thus enable high con- 
trast on unfriendly apertures . We first present a new approach 60 
to coronagraph design according to an embodiment of the 
current invention in the presence of a central obscuration, but 
in the absence of spiders or segments. We show that for 
coronagraphs with a pupil apodization and an opaque focal 
plane stop, contrasts of 10 -10 can be reached for any central 65 
obscuration diameter, provided that the Inner Working Angle 
is laige enough. Naturally the secondary support structures, 
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and in the segmented cases, segment gaps, will degrade this 
contrast. As our goal is to use two DMs as an amplitude 
modulation device, we first briefly review the physics of such 
a modulation. We then introduce a solution to tlii s problem. In 
particular, we show how to compute DM surfaces that miti- 
gate spiders and segment gaps. Current algorithms used for 
amplitude control operate under the assumption that ampli- 
mde errors are small, and thus they cannot be readily applied 
to the problem of compensating aperture discontinuities, 
which have inherently large reflectivity non-uniformities. 
FIG. 1G illustrates how the current example embodiment 
introduces a control strategy for the DMs that is radically 
different from previously published amplitude modulators in 
the high -contrast imaging literature. Ourtechnique, which we 
name Active Compensation of Aperture Discontinuities 
(hereafter ACAD), finds the adequate DM shapes in the true 
non-linear large amplitude error regime. In this case the DMs ’ 
surfaces are calculated as the solution of a non-linear partial 
differential equation, called the Monge-Ampere Equation. 
Below, we further describe our methodology according to an 
embodiment of the current invention to solve this equation 
and illustrate each step using an obscured and segmented 
geometry similar to JWST. As ACAD DM surfaces are pre- 
scribed in the ray-optics approximation, this is a fundamen- 
tally broadband technique provided that chromatic diffractive 
artifacts, edge ringing in particular, do not significantly 
impact the contrast. We find that when remapping small dis- 
continuities with Deformable Mirrors, the spectral bandwidth 
is only limited by wavelength-dependent edge-diffraction 
ringing in the Fresnel approximation (as discussed in Pueyo 
& Kasdin (2007) for instance). Below, we also present the 
application of some embodiments of the current invention to 
various observatory architectures. Some embodiments of the 
current invention can render high contrast coronagraphy pos- 
sible on any observatory geometry without adding any new 
hardware. 

2 Coronagraphy with a Central Obscuration 

2.1 Optimizing Pupil Apodization in the Presence of a Central 
Obscuration 

Because the pupil obscuration in an on-axis telescope is 
large it will be very difficult to mitigate its impact with DMs 
with a limited stroke. Indeed, the main hindrance to high- 
contrast coronagraphy in on-axis telescope is the presence of 
the central obscuration: it often shadows much more than 
10% of the aperture width while secondary supports and 
segments gaps cover -1%. We thus first focus of azimuthally 
symmetric coronagraphic designs in the presence of a central 
obscuration. This problem (without the support structures) 
has been addressed in previous publications either using cir- 
cularly symmetric pupil apodization (Sounnner et al. 201 1 ) or 
a series of phase masks (Mawet et al. 2011). Both solutions 
however are subject to limitations. The singularity at the 
center of the Optical Vector Vortex Coronagraph might be 
difficult to manufacture and a circular opaque spot thus lies in 
the central portion of the phase mask (<7/D) which results in 
a degradation of the ideal contrast of such a coronagraph 
(Krist et al. 2011). The solutions in Sounnner et al. (2011) 
result from an optimization seeking to maximize the off-axis 
throughput for a given focal plane stop diameter: the final 
contrast is absent from the optimization metric and is only a 
byproduct of the chosen geometry. Higher contrasts are then 
obtained by increasing the size of the focal plane mask, and 
thus result in a loss in IWA. 

2.2 The Optimization Problem 

Here we revisit the solution proposed by Sounnner et-al. 
(20 1 1 ) in a slightly different framework. We recognize that, in 
the presence of wavefront errors, high contrast can only be 
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achieved in an area of the field of view that is bounded by the 
spatial frequency corresponding to the DM’s actuator spac- 
ing. We thus consider the design of an Apodized Pupil Lyot 
Coronagraph which only aims at generating high contrast 
between the Inner Working Angle (IWA) and Outer Working 5 
Angle (OWA). In order to do so, we rewrite coronagraphs 
described by Soummer et~al. (2003) as an operator C which 
relates the entrance pupil P(r) to the electrical field in the final 
image plane. We first call P(lj) the Hankel transform of the 
entrance pupil: to 


0 ) 

where D is the pupil diameter, D s the diameter of the second- 
ary and § the coordinate at the science detector expressed in 
units of angular resolution (Xq/D). X 0 is the design wavelength 1 5 
of the coronagraph chosen to translate the actual physical size 
of the focal plane mask in units of angular resolution (often at 
the center of the bandwidth of interest). X is then the wave- 
length at which the coronagraph is operating (e.g. the physi- 
cal size of the focal plane mask remains constant as the width 20 
of the diffraction pattern changes with wavelength). For the 
purpose of the monochromatic designs presented herein 
X=X 0 . Then the operator is given by: 

25 

A 2 (2) 

C[«r)](f) = m - -j A; mt, n)ndn 

A 0 JD s /2 

where K(§,r|) is the convolution kernel that captures the effect 30 
of the focal plane stop of diameter M stop : 


m^)=[ 0 M ^ 2 J 0 (u^)J 0 (u%)udu 


(3) 


An analytical closed form for this kernel can be calculated 
using Lommel functions. Note that this Eq. 2 assumes that the 35 
Lyot stop is not undersized. Since we are interested in high 
contrast regions that only span radially all the way up to a 
finite OWA, we seek pupil apodization of the form: 


P(r) = 



40 


(4) 


where Ig(r) denotes the Bessel function of the first kind of 
order Q and a k e the kth zero of this Bessel function. In order 
to devise optimal apodizations over obscured pupils, Q, can 
be chosen to be large enough so that I e (r)«l for r<D^/2 (in 
practice we choose Q=1 0), The c% e corresponds to the spatial 50 
scale of oscillations in the coronagraph entrance pupil, and 
such a basis set yields high contrast regions all the way to 
OWA^N^^^Xq/D. Since the operator in Eq. 2 is linear, 
finding the optimal can be written as the following linear 
programming problem: 55 


max |min(P(r)) 1 Under the constraints 

(Pk) 1 r J 

(5) 

|C[(/’(r))l(f)| < 10-^ for IWA < f < OWA 

(6) 

ma x(P(r)) = 1 

(7) 

1 d 1 

(8) 

| — [f>(r)]| <b. 


Our choice of cost function and constraints has been 
directed by the following rationale: 

We maximize the smallest value on the apodization func- 
tion in an attempt to maximize throughput. The actual 
throughput is a quadratic function of the p k . Maximizing 
it requires the solution of a non-linear optimization 
problem (as described in Vanderbei et~al. (2003). 

The contrast constraint is enforced between the IWA and 
the OWA (<N modes ). 

The maximum of the apodization function is set to one 
(otherwise the p k will be chosen to be sufficiently small 
so that the contrast constraint is met). 

The absolute value of the derivative across the pupil cannot 
be larger than a limit, denoted as b here. As the natural 
solutions of such problems are very oscillatory (or “bang 
bang”, [Vanderbei et~al, (2003)], a smoothness con- 
straint has to be enforced (see [Vanderbei et~al. (2007)] 
for a similar case). 

Note that the linear transfer function in Eq. 2 can also be 
derived for other coronagraphs, with grayscale and phase 
image-plane masks, or for the case of under-sized Lyot stops. 
As general coronagraphic design in obscured circular geom- 
etries is not our main purpose, we limit the scope of the paper 
to coronagraphs represented by Eq. 2. 

Optimal design Apodized Pupil Lyot Coronagraphs on cir- 
cularly obscured apertures. With fixed obscuration ratio, size 
of opaque focal plane mask, IWA and OWA, our linear pro- 
gramming approach yields solutions with theoretical contrast 
below 1 0 10 . All PSFs shown on this figure are monochromatic 
for X=X 0 . When all other quantities remain equal and the 
central obscuration ratio increases (from 1 0% in the top panel 
to 20% in the middle panel), then the solution becomes more 
oscillatory (e.g less feasible) and the contrast constraint has to 
be relaxed. Eventually the optimizer does not find a solution 
and the size of the opaque focal plane mask (and thus the 
IWA) has to be increased (central obscuration of 30% in the 
bottom panel). On the right hand side, we present our results 
in two configurations: when the apodization is achieved using 
at grayscale screen (APLC), “on-sky” X/D bottom x-axis, and 
when the apodization is achieved via two pupil remapping 
mirrors (PIAAC), “on-sky” X/D top x-axis. We adopt this 
presentation to show that ACAD is “coronagraph indepen- 
dent” and that it can be applied to coronagraphs with high 
throughput and small IWA. 

2.3 Results of the Optimization 

Typical results of the monochromatic optimization in Eqs. 
8, with X=X 0 , are shown in FIG. 2 .2 for central obscurations of 
10, 20 and 30%. In the first two cases the size of the focal 
plane stop is equal to 3X/D, the IWA is 4X/D and the OWA is 
30X/D. As the size of the central obscuration increases the 
resulting optimal apodization becomes more oscillatory and 
the contrast constraint has to be loosened in order for the 
linear programming optimizer to converge to a smoother 
solution. Alternatively increasing the size of the focal plane 
stop yields smooth apodizers with high contrast, at a cost in 
angular resolution (bottom panel with a central obscuration of 
30%, a focal plane mask of radius 4X/D, an IWA of 5X/D and 
an OWA of 30X/D). These trade-offs were described in Soum- 
mer et~al. (201 1 ), however our linear programming approach 
to the design of pupil apodizations now imposes the final 
contrast instead of having it be a by product of fixed central 
obscuration and focal plane stop. These apodizations can 
either be generated using grayscale screen (at a cost in 
throughput and angular resolution) or a series of two aspheri- 
cal PIAA mirrors (for better throughput and angular resolu- 
tion). In order not to lose generality, we present our results on 
FIG. 2.2 considering the two types of practical implementa- 
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tions (classical apodization and PIAA apodization). In the 
case of a grayscale amplinide screen the angular resolution 
units are as defined in Eq. 2 and the throughput is smaller than 
unity. In the case of PIAA apodisation the throughput is unity 
and the angular resolution units have been magnified by the 
field independent centroid based angular magnification 
defined in Pueyo et~al. (201 1 ). We adopt this presentation for 
the remainder of the paper where one dimensional PSFs will 
be presented with “APLC angular resolution units” in the 
bottom horizontal axis and“PIAAC angular resolution units” 
in the top horizontal axis. 

Note that this linear programming approach only optimizes 
the contrast for a given wavelength. However, since the solu- 
tions presented in FIG. 2.2 feature contrasts below 10 10 , we 
choose not to focus on coronagraph chromatic optimizations. 
Instead, in order to account for the chromatic behavior of the 
coronagraph, the monochromatic simulations in §. 6 are car- 
ried out under the conservative assumption that the physical 
size of the focal plane stop is somewhat smaller than optimal 
(or that the operating wavelength of the coronagraph is 
slightly off, X=1.2X 0 ). As a consequence the raw contrast of 
the coronagraphs presented in §. 6 is ~10~ 9 . Note that this 
choice is not representative of all possible Apodized Pupil 
Coronagraph chromatic configurations . It is merely a shortcut 
we use to cover the variety of cases presented in §. 6. In §. 7.2 
we present a set of broadband simulations that include wave- 
front errors and the true coronagraphic chromaticity for a 
specific configuration and show that bandwidth is more likely 
to be limited by the spectral bandwidth of the wavefront 
control system than by the coronagraph. However, future 
studies aimed at defining the true contrast limits of a given 
telescope geometry will have to rely on solutions of the linear 
problem in Eqs. 8 which has been augmented to accommo- 
date for broadband observations. In theory, the method pre- 
sented here can also be applied to asymmetric pupils. How- 
ever, the optimization quickly becomes computationally 
intensive as the dimensionality of the linear programming 
increases (in particular when the smoothness constraint and 
the bounds on the apodization have to be enforced at all points 
of a two dimensional array). This problem can be somewhat 
mitigated when seeking for binary apodizations, as shown in 
Carlotti et~al. (2011), at a cost in throughput and angular 
resolution. 

3 Physics of Amplitude Modulation 
3.1 General Equations 

We have shown in §. 2 that by considering the design of 
pupil apodized coronagraphs in the presence of a circular 
central obscuration as a linear optimization problem, high 
contrast can be reached provided that the focal plane mask is 
large enough. In practice, the secondary support structures 
and the other asymmetric discontinuities in the telescope 
aperture (such as segment gaps) will prevent such levels of 
starlight suppression. We demonstrate that well controlled 
DMs can circumvent the obstacle of spiders and segment 
gaps. In this section, we first set-up our notations and review 
the physics of phase to amplitude modulation. We consider 
the system represented on FIG. 3 where two sequential DMs 
are located between the telescope aperture and the entrance 
pupil of the coronagraph. In this configuration, the telescope 
aperture and the pupil apodizer are not in conjugate planes. 
This will have an impact on the chromaticity of the system 
and is discussed in §. 5. Without loss of generality we work 
underthe “folded” assumption illustrated on FIG. 3 where the 
DMs are not tilted with respect to the optical axis and can be 
considered as lenses of index of refraction -1 (as discussed in 
[Vanderbei Traub(2005)]. In the scalar approximation the 


relationship between the incoming field, E DM1 (x,y), and the 
outgoing field, E Z3iU2 (x 2 ,y 2 ), is given by the diffractive Huy- 
gens Integral: 


5 

1 f .2ir„, , (9) 

Edm 2 (x 2 , y2) = T-g j E dm , (.*, y)e X e<W2J ' 2) dxdy 

where <A- corresponds to the telescope aperture and Q(x,y, 
x 2 ,y 2 ) stands for the optical path length between any two 
points at DM1 and DM2: 

Q(x,y,x 2 ,y 2 )'=h i (x,y)+S(x,y,X2,y2)-h2(X2.y2) ( 10 ) 

S(x,y,x 2 ,y 2 ) is the free space propagation between the DMs: 




-x 2 f + (.v -y 2 ) 2 + (Z + y) - li 2 {x 2 , y 2 )) 2 


( 11 ) 


where Z is the distance between the two DMs, tq and h 2 are 
the shapes of DM1 and DM2 respectively (as shown on FIG. 
3) and k is the wavelength. We recognize that two sequential 
DMs act as a pupil remapping unit similar to PIAA corona- 
graph [Guyon(2003)] whose ray optics equations were first 
derived by [Traub & Vanderbei (2003)]. We briefly state the 
notation used to describe such an optical system as introduced 
in Pueyo et~al. (2011): 

For a given location at DM2, (x 2 ,y 2 ), the location at DM1 
30 of the incident ray according to ray optics is given by: 

x 1 (x 2 ,y 2 )=f l (x 2 ,y 2 ) (12) 

yi(x 2 ,y 2 )=Si(x 2 ,y 2 )- (13) 

Conversely, for a given location at DM1, (x^yj, the loca- 
tion at DM2 of the outgoing ray according to geometric 
optics is given by: 

x 2 (x l ,y I )=f 2 (x„y 1 ) (14) 

40 y2(xi,y,)=g 2 (x 1 ,y 1 ). ( 15 ) 

Fermat’s principle dictates the following relationships 
between the remapping functions and the shape of DM1 : 


dh { 

- fL 

-f 2 (x l ,y i ) 

(16) 

dx 


z 


dh { 

- — 

-gl&i, yi) 

(17) 

dy 

lUlD'l) 

z 



Conversely, if we choose the surface of DM2 to ensure that 
the outgoing on-axis wavefront is flat, we then find: 


dh 2 


x 2 -fi(x 2 ,y 2 ) 

(18) 

dx 

'(x 2 ,y 2 ) 

z 


dh 2 

| 

y2 - gi (*2, y2) 

(19) 

dy 

'(x 2 ,y 2 ) 

z 



Schematic of the optical system considered: the telescope 
apertures is followed by two sequential Defromable Mirrors 
(DMs) in non-conjugate planes whose purpose is to remap the 
65 pupil discontinuities. The beam then enters a coronagraph to 
suppress the bulk of the starlight: in this figure we show an 
Apodized Lyot Pupil Coronagraph (APLC). This is the coro- 
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nagraphic architecture we consider for the remainder of the 
paper but we stress that the method presented herein is appli- 
cable to any coronagraph. 

3.2 Fresnel Approximation and Talbot Imaging 

In Pueyo et~al. (201 1) we showed that one could approxi- 5 
mate the propagation integral in Eq. 9 by taking in a second 
order Taylor expansion of Q(x,y,x 2 ,y 2 ) around the rays that 
trace (fi(x 2 ,y 2 ), gi(x 2 ,y 2 )) to (x 2 ,y 2 ). In this case the relation- 
ship between the fields at DM1 and DM2 is: 

to 


EpM 2 (*2* >2) = 


(20) 


e'T 

iXZ 


\L 


E D m i (ac, y)e 


• 2+2 'a^ 0-^1 
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dxdy \ 




Wlien the mirror’s deformations are very small compared 20 
to both the wavelength and D 2 /Z, the net effect of the wave- 
front disturbance created by DM1 can be captured in E Z)M1 
(x,y) and the surface of DM2 can be factored out of Eq. 9. In 
this case 


*i =x 2 , y i =y 2 , 


dh 

dx 


= 1 , 


ah 

dy 


= 0 . 


dg 2 
dy 


= 1. 


Then, Eq. 20 reduces to: 


30 


£dM2(*2> >2) = 

oz~ 


( 21 ) 


f c x'*i Mexzfr-iiWy-yz? ) dxdy 

J Aperture 
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structures and segments) with the DMs. In practice, when 
using Eq. 22 in the wavefront control scheme outlined in 
Pueyo et~al. (2009) to correct aperture discontinuities, this 
weak coupling results in large mirror shapes that lie beyond 
the range of the linear assumption made by the DM control 
algorithm. For this reason, methods outlined on the left panel 
of FIG. 1G to correct for aperture discontinuities do not 
converge to high contrast. Because phase to amplitude con- 
version is fundamentally a very non-linear phenomena, these 
descending gradient methods [Borda & Traub (2006), 
Give’ on et-al. (2007) Pueyo et~al. (2009)] are not suitable to 
find DM shapes that mitigate apertures discontinuities. We 
circumvent these numerical limitations by calculating DMs 
shapes that are based on the full non-linear problem, right 
panel of FIG. 1G. 

3.3 The SR-Fresnel Approximation 

In the general case, starting from Eq. 20 and following the 
derivation described in §. 5, the field at DM2 can be written as 
follows: 


£ d « 2 (.* 2 , yi) = | V|det[7]| 
Jtp 


(23) 


K*2>?2) 


where S DMi (S.,r| ) is the Fourier transform of the telescope 
aperture, T T and stands for the Fourier plane. We call this 
integral the Stretched-Remapped Fresnel approximation 
(SR-Fresnel). Moreover det[J(x 2 ,y 2 )] is the determinant of 
the Jacobian of the change of variables that maps (x 2 ,y 2 ) to 


which is the Fresnel approximation. If moreover 


40 


det[7U 2 , 3*2 )] = 


dfi dg, 
dx dy 



(24) 


(2 7t \ 

«i (■*> y) = A^cosl — (mx + ny)\, 
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h 2 (x,y)=-hj (x,y), with e«l , then the outgoing field is to first 
order: 


7 tA Z( m 2 +n 2 ) j In A (22) 50 

Edmi(x 2 , y 2 ) — 2 .IecosI — (»i.v 2 + ny 2 ) I. 

nAZ(m 2 +n 2 ) (2ji A 

Edmi(x 2 , y 2 ) = — 2 Aecos| — (»i.v 2 +ny 2 ) j 
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Tliis phase-to-amplitude coupling is a well known optical 
phenomenon called Talbot imaging and was introduced to the 
context of high contrast imaging by [Shaklan & Green 
(2006)]. In the small deformation regime, the phase on DM1 
becomes an amplitude at DM2 according to the coupling in 60 
Eq. 22. Wlien two sequential DMs are controlled to cancel 
small amplitude errors, as in Pueyo et-al. (2009), they operate 
in this regime. Note, however, that the coupling factor scales 
with wavelength (the resulting amplitude modulation is 
wavelength independent, but the coupling scales as X): this 65 
formalism is thus not applicable to our case, for which we are 
seeking to correct large amplitude errors (secondary support 


In the ray optics approximation, X-0, the non-linear trans- 
fer function between the two DMs becomes: 


Eu^/ X j..yj) { ^\del/./\i\ (25) 

\E-uhn ( x 2 jA 1 2 = { det[J] [E DMl (26) 

The square form (e.g Eq. 26) of this transfer function can 
also be derived based on conservation of energy principles 
and is a generalization to arbitrary geometries of the equation 
driving the design of PIAA coronagraphs [Vanderbei & Traub 
(2005)]. A full diffractive optimization of the DM surfaces 
requires use of the complete transfer function shown in Eq. 
23. However, there do not exist yet tractable numerical 
method to evaluate Eq. 23 efficiently enough in order for tliis 
model to be included in an optimization algorithm. Moreover 
even solving the ray optics problem is extremely compli- 
cated: it requires to find the mapping function (f, ,gj ) which 
solves the non-linear partial differential equation in Eq. 26. 
Substituting for (fj ,gj ) and using Eqs. 19 yields a second 
order non-linear partial differential equation in h 2 . This is the 
problem that we set ourselves to tackle in the next section, and 
is the cornerstone of our Adaptive Compensation of Aperture 
Discontinuities. As a check, one can verify that in the small 
deformation regime 



23 


US 9,140,896 B2 


24 


e.g. if hi (. x , y ) = 


. - j — (mx + ny ) j and h 2 (x, y) = -hi(x, y)j 


[ J df-, dg-> 



[G)Ml(*l> y ! )] 2 


(30) 


Eq. 26 yields the same phase-to-amplitude coupling as in 
Talbot imaging [Pueyo (2008)], Eq. 26 is a well know optimal 
transport problem [Monge (1781)], which has already been 
identified as underlying optical illumination optimizations 
[Glimm & Oliker (2002)]. While the existence and unique- 10 
ness of solutions in arbitrary dimensions have been exten- 
sively discussed in the mathematical literature (see [Dacoro- 
gna & Moser (1990)] for a review), there was no practical 
numerical solution published up until recently. In particular, 
to our knowledge, not even a dimensional solution for which 15 
the DM surfaces can be described using a realistic basis-set 
has been published yet. We now introduce a method that 
calculates solutions to Eq. 26 which can be represented by 
feasible DM shapes. 

4 Calculation of the Deformable Mirror Shapes 20 

4.1 Statement of the Problem 

Ideally, we seek DM shapes that fully cancel all the dis- 
continuities at the surface of the primary mirror and yield a 
uniform amplitude distribution, as shown in the top panel of 
FIG. 2. A solutions for a particular geometry with four sec- 25 
ondary support has been derived by Lozi et~al. (2009). It 
relies on reducing the dimensionality of the problem to the 
direction orthogonal to the spiders. It is implemented using a 
transmissive correcting plate that is a four-faced prism 
arranged such that the vertices coincides with the location of 30 
the spiders. The curvature discontinuities at the location of the 
spiders are responsible for the local remapping that removes 
the spiders in the coronagraph pupil. However such a solution 
cannot be readily generalized to the case of more complex 
apertures, where the secondary support structures might vary 35 
in width, or in the presence of segment gaps. Moreover it is 
transmissive and thus highly chromatic. Here we focus on a 
different class of solutions and seek to answer a different 
question. How well can we mitigate the effect of pupil dis- 
continuities using DMs with smooth surfaces, a limited num- 40 
her of actuators (e.g a limited maximal curvature), and a 
limited stroke? Under these constraints directly solving Eq. 

26 (e.g. solving the forward problem illustrated in the top 
panel of FIG. 2) is not tractable as both factors on the left hand 
side of Eq. 26 depend on h 2 . More specifically, the implicit 45 
dependence of ^. DMl (fi(x 2 ,y 2 )- gife^)) on h 2 can only be 
addressed using finite elements solvers, whose solutions 
mi gilt not be realistically representable using a DM. However 
this can be circumvented using the reversibility of light and 
solving the reverse problem, where the two mirrors have been 50 
swapped. Indeed, since 

X2=f2ifl(X2,y 2 )Sl(x 2 .y2)) ( 27 ) 


Since our goal is to obtain a pupil as uniform a possible we 
seek a field at DM2 as uniform as possible: 

Gm/>( 

V a tE I>Ml (x,y) 2 dxdy=Constmt (31) 

Moreover we are only interested in compensating asym- 
metric structures located between the secondary and the edge 
of the primary. We thus only seek to find (f 2 ,g 2 ) such that: 


Edm i(*i, y i) = 


(32) 


Mxi, .vi) = l[P(.v. y) - (1 - P c (x, y))]*e 




■(*101) 


where P 0 (x 1 ,y 1 ) is the obscured pupil, without segments or 
secondary supports. Finally, we focus on solutions with a 
high contrast only up to a finite OWA. We artificially taper the 
discontinuities by convolving the control term in the Monge- 
Ampere Equation, [P(x,y)-(1-P 0 (x,y)], with a gaussian of 
width co, Note that this tapering is only applied when calcu- 
lating the DM shapes via solving the reverse problem. When 
the resulting solutions are propagated through Eq. 26 we use 
the tme telescope pupil for E zwl (x 1 ,y 1 ). The parameter to has 
a significant impact on the final post-coronagraphic contrast. 
Indeed we are here working with a merit function that is based 
on a pupil-plane residual, while ideally our cost function 
should be based on image-plane intensity. By convolving the 
control term in the reverse Monge- Ampere Equation, we low 
pass filter the discontinuities. This is equivalent to giving a 
stronger weight to low-to-mid spatial frequencies of interest 
in the context of exo-planets imaging. For each case presented 
in §. 6 we calculate our DM shapes over a grid of values of to 
which correspond to low pass filters with cutoff frequencies 
ranging from -OWA to -2 OWA and we keep the shapes 
which yield the best contrast. This somewhat ad-hoc 
approach can certainly be optimized for higher contrasts. 
However such an optimization is beyond the scope of the 
present manuscript. 

The problem we are seeking to solve is illustrated in the 
second panel of FIG. 2. In this configuration the full second 
order Monge- Ampere Equation can be written as: 



1+z 


a 2 hi 

' dy 2 



= Mx, y ) 2 


(33) 


55 


y2=g2(fi(X2,y2)&l(X2,y2)) (28) 

then we have the following relationship between the determi- 
nants of the forward and reverse remappings: 

. _ ISfz dgi (%U : 1 I ( d/i dgi tdgi T] I 
~\dx dy Vdx) j \ dx dy { dx ) j l^, 

We thus focus on the inverse problem, bottom panel of FIG. 65 
2, that consists of first finding the surface oTh , as the solution 
of: 


where we have dropped the (x, ,y , ) dependence for clarity. 
Since we are interested in surface deformations which can 
realistically be created using a DM, we seek for a Fourier 
representation of the DMs surface: 

D 2 (34) 

hi(x,y)=—HUX, Y) 

i,2 N/2 m 

= - Z Z a ^ (mX+nX) 

Z n=-N/2m=-N/2 


60 

(29) 
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with 


a -m,-n = a -m-n > a -m,-n = Q-m-n, 

where N is the limited number of actuators across the DM. 
Note that we have normalized the dimensions in the pupil 
plane X=x/D, Y=y/D. The nonnalized second order Monge- 
Ampere Equation is then: 


( 1 + 



d 2 Hi 1 ( 

dY 2 ) \3XdYl 


= A(X, Yf 


(35) 


15 


For each configuration in this paper we first solve Eq. 35 
and then transform the normalized solution in physical units, 
which depends on the DMs diameter D and their separation Z. 

4.2 Solving the Monge- Ampere Equation to Find H t 20 

Over the past few years a number of numerical algorithms 
aimed at solving Eq. 35 have emerged in the literature 
[Loeper & Rapetti (2005), Benamou et~al. (2010)]. Flere we 
summarize our implementation of two of them: an explicit 
Newton method [Loeper &Rapetti (2005)], and a semi-im- 25 
plicit method [Froese & Oberman (2012)]. We do not delve 
into the proof of convergence of each method, they can be 
found in Loeper & Rapetti (2005), Benamou et~al. (2010) 
Froese & Obennan (2012)]. Note that Zheligovsky et~al. 3(J 
( 2010 ) discussed both approaches in a cosmological context 
and devised Fourier based solutions. Flere we are interested in 
a two dimensional problem and we outline below the essence 
of each algorithm. 

4.2.1 Explicit Newton Algorithm 35 

This method was first introduced by [Loeper & Rapetti 
(2005)] and relies on the fact that Eq. 35 can be re-written as 


[( d 2 Hi 

1 + IF 

d 2 Hi 

dXdY 

d 2 H, 

, a 2 Hi 

1+ BY 2 

LI dXdY 


= det [ld+'H[Hi(X, F))] 


(36) 40 


L E (H k ,V k ) = (39) 

d 2 H k l )d 2 V t ( t d 2 Hi)d 2 V k „d 2 //[ d 2 V k 
{ + dY 2 jdX 2 + dX 2 J dY 2 ~~dXdY dXdY 

, 1 , . (40) 

Re(H\)= -(A 2 — det[W+'W(//[)]) 


and solve 

Le{.H*, V k )=R E (Hi k ). (41) 

Eq. 41 is a linear partial differential equation in V*. Since 
we are interested in a solution which can be expanded in a 
Fourier series we write V k as: 


m nh (42) 

Vf(X, Y) = jr jr v k m .„e’T>' mX+ ’ ,X) . 

n=-NI 2 m=—Nj 2 


Both the riglat hand side and the left hand side of Eq. 45 can 
be written as a Fourier series, with a spatial frequency content 
between -N and N cycles per aperture. Equating each Fourier 
coefficient in these two series yields the following linear 
system of (2N+1 ) 2 equations with (N+l ) 2 unknowns. 


For all m 0 ,n 0 t[-N,N]:j DM e ankm,>x +" ,>r> [L E (Hi k , t*)- 

R E (H l k )]dXdY= 0 (43) 

When searching for V* as a Fourier series over the square 
geometry chosen here, this inverse problem is always well 
posed. 

Update the solution 

The convergence of this algorithm relies on the introduc- 
tion of a damping constant x>l . [Loeper & Rapetti (2005)] 
showed that as long as 


X 2 + r 


- + wf 


where W (.) is the two dimensional Hessian of a scalar field 
and Id the identity matrix. If one writes H 3 =u+v with ||v||«||u|| 
then: 

det[Id+ Tt ( u+&v)]=det[Id+H(ii)]+bTr[{Id+H(ii)f T H 

Ml+olS 2 ) (37) 

where (,)' r denotes the transpose of the comatrix. Eq. 35 can 
thus be linearized as: 

( d 2 u\d 2 v ( d 2 u\d 2 v n d 2 u d 2 v (38) 

( + 3F 2 JaX 2+ [ + dX 2 jdY 2 ~~dXdY dXdY = 

^(X, Y) 2 - det^U-^L + 

Hie explicit Newton algorithm relies on Eq. 38 and can 
then be summarized as carrying out the following iterations: 
Choose a first guess Hj°. 

At each iteration k we seek for a solution of the form 
IU + ' IU+V' 1 , where W k is the DM shape update. 

In order to find V k we write: 


remains convex, which is always true for ACAD with reason- 
ably small aperture discontinuities, there exists a x large 
enough so that this algorithm converge towards a solution of 
Eq. 35. However since this algorithm is gradient based, it is 
not guaranteed that it converges to the global minimum of the 
underlying non-linear problem. In order to avoid having this 
solver stall in a local minimum we follow the methodology 
outlined by [Froese & Oberman (2012)] and first carry out a 
series of implicit iterations to get within a reasonable neigh- 
borhood of the global minimum. 

Virtual field at Ml when solving the reverse problem. A 
(XjjYj) is the desired apodization (top), A"(X 1 ,Y 1 ) is the 
apodization obtained after solving the Monge-Ampere Equa- 
tion (center). The bottom panel shows the difference between 
the two quantities: the bulk of the energy in the residual is 
located in high spatial frequencies that cannot be controlled 
by the DMs. 

4.2.2 Implicit Algorithm 

This algorithm, along with its convergence proof, is thor- 
oughly explained in [Froese Oberman (2012)]. It relies on 
rewriting Eq. 35 as: 


45 


50 
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a 2 ^ a 2 Hi I ; — (44) 

~av + = v det[,d + ww > (X - D)l 2 +24,A '’ Y) 

The implicit method consists of carrying out the following 
iterations: 

Choose a first guess H, 0 
In order to find FU 1 we write: 


V k ) = 


a 2 H k i +l 

ax 2 


a 2 H k i +1 

av 2 


R E (H k ) = -y/det [ld + WHi(X, Y))f + 2A(X, Y ) 2 


and solve 

L 1 (H l M )=R l (Hi k ). (45) 

This problem is a linear system of (N+l) 2 equations with 
(N+l) 2 unknowns and can be solved using projections on a 
Fourier Basis: 

For all m o ,n oE [-N,N]:S DM e a ’' (m0X *''’ ,Y> [Li(Hi k , J*)-rt, 

(Hf)]dXdY=0 (46) 

Note that the term under the square root in is 

guaranteed to be positive at each iteration. 

Iterate over k 

The inverse problem in Eq. 46 is always well posed, for any 
basis set or pupil geometry, while the explicit Newton method 
runs into convergence issues when not using a Fourier basis 
over a square. When seeking to use a basis set that is more 
adapted to the geometry of the spiders and segments or when 
using a trial influence function basis for the DM, the implicit 
method is the most promising method. In this paper we have 
limited our scope to solving the reverse problem in the bottom 
panel of FIG. 4 , and used a Fourier representation for the DM, 
we are able to use both methods. In order to make sure that the 
algorithm converges towards the true solution of Eq. 35 we 
first rim a few tens of iterations of the implicit method and, 
once it has converged, we seek for a more accurate solutions 
using the Newton algorithm. Typical results are shown in 
FIG. 6 where most of the residual error resides in the high 
spatial frequency content (e.g. above N cycles per aperture). 
Our solutions are limited by the non-optimality of the Fourier 
basis to describe the mostly radial and azimuthal structures 
present in telescopes’ apertures. Moreover the DM shape is 
the result of the minimization of a least squares residual in the 
virtual end-plane of the reverse problem, with little regard to 
the spatial frequency content of the solution in the final image 
plane of the coronagraph. While this method yields signifi- 
cant contrast improvements, as reported in §. 6, we discuss in 
§. 7 how it can be refined for higher contrast. 

4.3 Deformation of the Second Mirror 

Once the surface of DM1 has been calculated as a solution 
of Eq. 26, we compute the surface of DM2 based on Eqs. 19, 
which stem from enforcing flatness of the outgoing on-axis 
wavefront. We seek a F ourier representation for the surface of 
DM2: 


F 2 (X l ,Y l )=X l - 


£ i2n(mXi+nYi) 


G 2 (X h Yi)=Yi- 

2 M 2J. M2--,.- n J2x(f*Xn-nYi) 

5 ■ -M2 a mji e 

Then the normalized version of Eqs. 19 can be rewritten as: 
L x (Xi,Yi)= 

10 

R x (Xi,Yi)-Xi-F 2 (Xi,Yi) 

L X =R X (48) 

15 L y (Xi,Yi)= 

j:^_ Na Na i2mb m ^^ x ‘- y ‘^ kx ‘- r m 

R y (X l ,r l )=Y l -G 2 (X l ,Y l ) 

20 L=R y (49) 

We then multiply each side of Eq. 48 and Eq. 49 by: 

e i2xim^2aW+noF2iXi.Yl)) det p d+ ft (pk^Y^)] 

where (m 0 ,n 0 ) corresponds to a given DM spatial frequency. 
22 Integrating over the square area of the DM and using the 
orthogonality of the Fourier basis yields the following system 
of 2* (N+ 1 ) 2 equations with (N+ 1 ) 2 real unknowns : 


30 


35 


40 


For all(/«o, «o): 

2 m m 0 b m()/l0 = 

I R x (X. Y)det[ld + 'H(H l i(X, Y))]e BM ™> F2kx - r * n <> F * x ™dXdY 
Jdm 

For all(/«o, no): 

Uni n 0 bm 0 M 0 = 

f R x (X, Y)Aet[ld + 'H(H k i(X, y))] c i2r < m o F 2<*.r)+"o F 2<*4)) dM >< 
Jdm 

2ninob mo ,„ 0 = 

f R y (X, Y)dtt[ld + mH k i(X , Y))]e‘ 23t(m <> F2(x ’ r)+ ' , o F2(x ’ r)) dXdY 
Jdm 


45 

We then find H 2 , the normalized surface of DM2, by solving 
this system in the least squares sense. Once the Monge- 
Ampere Equation has been solved, the calculation of the 
surface of the second mirror is a much easier problem. Indeed, 
50 by virtue of the conservation of the on-axis optical path 
length, finding the surface of DM2 only consists of solving a 
linear system (see [Traub & Vanderbei (2003)]). 

4.4 Boundary Conditions 

The method described above does not enforce any bound- 
55 ary conditions associated with Eq. 35. One practical set of 
boundary conditions consists of forcing the edges of each DM 
to map to each other: 


hi(x, y ) = —H 2 (X,Y) 


42 2* 

n=—N/2 m=-N/2 


^-^■(/nX+nX) 


(47) 


Plugging the solution found in the previous step for h t into 65 
Eqs. 17 yields a closed form for the normalized remapping 
functions, (F 2 ,G 2 ): 


444 


( 50 ) 

( 51 ) 

( 52 ) 
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-continued 



with i=l,2. These correspond to a set of Neumann boundary 
conditions in Hj(X,Y) and H 2 (X,Y). These bomidary condi- 
tions can be enforced by augmenting the dimensionality of 
the linear systems on Eq. 43 and Eq. 46, however doing so 
increases the residual least squares errors and thus hampers to 
the contrast of the final solution. Moreover FIG. 4 and FIG. 5 
show that, because of the one to one remapping near the DM 
edges in the control term of the reverse problem, the boundary 
conditions are almost met in practice. For the remainder of 
this paper we thus do not include boundary conditions when 15 
calculating the DM shapes, when solving for Hj(X,Y) in Eq. 

35 since, in the worse case, only the edge rows and columns 
of the DMs actuator will have to be sacrificed in order for the 
edges to truly map to each other. 

4.5 Remapped Aperture 20 

For a given pupil geometry we have calculated (H, ,H 2 ). We 
then convert the DM surfaces to real units, (hj,h 2 ), by multi- 
plication with D 2 /Z. We evaluate the remapping functions 
using Eqs . 1 7 and 1 9 and obtain the field at the entrance of the 
coronagraph in the ray optic approximation 25 
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over large spectral windows. However, when taking into 
account the edge diffraction effects that are captured by the 
quadratic integral in Eq. 23, the true propagated field at DM2 
becomes wavelength dependent. More specifically, when k is 
not zero then the oscillatory integral superposes on the ray 
optics field a series of high spatial frequency oscillations. In 
theory, it would be best to use this as the full transfer function 
to include chromatic effects in the computation of the DMs 
shapes. However, as discussed in §. 4, solving the non-linear 
Monge-Ampere Equation is already a delicate exercise, and 
we thus have limited the scope of this paper to ray optics 
solutions. Nonetheless, once the DMs’ shapes have been 
determined using ray optics, one should check whether or not 
the oscillations due to edge diffraction will hamper the con- 
trast. This approach is reminiscent of the design of PIAA 
systems where the mirror shapes are calculated first using 
geometric optics and are then propagated through the diffrac- 
tive integral in order to check a posteriori whether or not the 
chromatic diffractive artifacts are below the design contrast 
[Pluzhnik et~al. (2005)]. In this section we detail the deriva- 
tion of Eq. 23 that is the diffractive integral for the two DMs 
remapping system and use this formulation to discuss the 
diffractive properties of ACAD. 

We start with the expression of the second order diffractive 
field at DM2 as derived in [Pueyo et~al. (2011)], Eq. 20. 


EdM2(X2, }'2) = 


(54) 


{VdSfTT 30 E DM2 (x 2 , y 2 ) = 


where the exponential factor corresponds to the Optical Path 
Length through the two DMs. Even if the surface of the DMs 
has been calculated using only a finite set of Fourier modes, 35 
we check that the optical path length is conserved. 

FIG. 6 shows that since the curvature of the DMs is limited 
by the number of modes N, our solution does not hilly map 
out the discontinuities induced by the secondary supports and 
the segments. However, they are significantly thinner and one 40 
can expect that their impact on contrast will be attenuated by 
orders of magnitude. In order to quantify the final corona- 
graphic contrasts of our solution we then propagate it through 
an APLC coronagraph designed using the method in §. 2. In 
the case of a hexagon based primary (such as JWST), we use 45 
a coronographic apodizer with a slightly oversized secondary 
obscuration and undersize outer edge in order to circularize 
the pupil. Note that this choice is mainly driven by the type of 
coronagraph we chose in § . 2 to illustrate our technique. Since 
the DM control strategy presented in this section is indepen- 50 
dent of the coronagraph, it can be generalized to any of the 
starlight suppression systems which have been discussed in 
the literature. For succinctness we present our results using 
coronagraph solely based on using pupil apodization (either 
in an APLC orinaPIAAC configuration). Results for a JWST 55 
geometry are shown on FIG. 5.3 and FIG. 7 and discussed in 
§. 6. Eq. 54 assumes that the propagation between the two 
DMs occurs according to the laws of ray optics. In the next 
section we derive the actual diffractive field at DM2, e.g Eq. 

23, and show that in the pupil remapping regime of ACAD, 60 
edge ringing due to the free space propagation is actually 
smaller than in the Fresnel regime. 

5 Chromatic Properties 

5.1 Analytical Expression of the Diffracted Field 

ACAD is based on ray optics. It is an inherently broadband 65 
technique, and provided that the coronagraph is optimized for 
broadband performance ACAD will provide high contrast 
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We write E DM1 (x,y) as its inverse Fourier transform 

E DMl (x,y)=l£ DMl ^)e‘ 2 ’’<^>c^dr\ (55) 

and insert this expression in Eq. 20. Completing the squares in 
the quadratic exponential factor then yields: 


Edm 2 U 2 , y 2 ) = --- 
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F/jMI (f. 


(56) 
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with 


h(yi,xi) = 



|A Z ,~Bg 2( v 1,^2 
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*dx 


(57) 



The integral over space, I , (y , ,x , ), is the integral of a com- 
plex gaussian and can be evaluated analytically. This yields: 
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EdMi(x 2 , V2) = 1 


/ df 2 dg 2 

fSg2 \ 2 

V dxdy 

\ <9jc 


f Edui (f , 


(myi) 


We thus have expressed E DA/2 (x 2 ,y 2 ) as a function of (x x , 
y 1 )=(f 1 (x 2 ,y 2 ),g 1 (x 2 ,y 2 )). This expression can be further sim- 
plified: using Eqs. 27 to 29 one can derive 


Bfi | 1 I 

dx W,!) - det[J(x 2 ,y 2 )] dy | (X2J , 2) 

and 

dg2\ 1 dfi I 

dy det[y(^2, y2)] d* I (X 2 ,y 2 ) 

Which finishes to prove Eq. 23: 
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EdM 2(*2, y2) = 

|vjd«yjr f 




(59) 
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l(*2 0>2 ) 


This expression is very similar to a modified Fresnel propa- 
gation and can be rewritten as such: 


35 


EdM2(x 2 , >2) = 


( 60 ) 
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Edmi{x 2, y2) = I i (*> y)e 
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Because of this similarity we call this integral the 
Stretched-Remapped Fresnel approximation (SR-Fresnel). 55 
Indeed in this approximation the propagation distance is 
stretched by 


a Jl a A 60 

By dx 

,det[y]’ det[J], 

and the integral is centered around the remapped pupil (f x ,gj ). 

5.2 Discussion 65 

FIG. 10 , Left: Comparison of edge diffraction between 
PIAA and Fresnel propagation. Because the discontinuities in 
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the pupil occur at the location where T V <1 , a PIAA amplifies 
the chromatic ringing when compared to a more classical 
propagation. Right: Comparison of edge diffraction between 
ACAD and Fresnel propagation. Because the discontinuities 
in the pupil occur at the location where T X >1 , a ACAD damps 
the chromatic ringing when compared to a more classical 
propagation. These simulations were carried out for a one 
dimensional aperture. In this case the Monge- Ampere Equa- 
tion can be solved using finite elements and the calculation of 
the diffractive effects reduces to the evaluation of various 
Fresnel special functions at varying wavelength. The full 
two-dimensional problem is not separable and requires the 
development of novel numerical tools. 

The integral form provides physical insight about the 
behavior of the chromatic edge oscillations. If we write 


r. 


det[J] 

dgi 

dy 


1**2 ->2 > 


_ det[7] 

d_h 

dx 


C2.V2) 


( 61 ) 


( 62 ) 


we can identify a several diffractive regimes: 

1. when T X =T =1 the E Z5AG (x,y 2 ) reduces to a simple 
Fresnel propagation. Edge ringing can them be miti- 
gated using classical techniques such as pre-apodiza- 
tion, or re-imaging into a conjugate plane using over- 
sized relay optics. 

2. when F A , T^cl then it is as if the effective propagation 
length through the remapping unit was increased. This 
magnifies the edge chromatic ringing. 

3. when I j. T >1 then it is as if the effective propagation 
length through the remapping unit was decreased. This 
damps the edge chromatic ringing. 

4. when r A >l and T r <l , e.g. at a saddle point in the DM 
surface, then it is as if the effective propagation length 
through the remapping unit was decreased in one direc- 
tion and increased in the other. The edge chromatic 
ringing can either be damped or magnified depending on 
the relative magnitude of T x and T . 

In the case of a PIAA coronagraphs, the mirror shapes are 
such that T x , r,.>l at the center of the pupil and r A .,r j ,«l at 
the edges of the pupil, where the discontinuities occur. As a 
consequence the edge oscillations are largely magnified when 
compared to Fresnel oscillations (see right panel of FIG. 10 ), 
and apodizing screens are necessary in order to reduce the 
local curvature of the mirror’s shape (as also discussed in 
[Pluzhnik et~al. (2005), Pueyo et~al, (201 1 )] . In the case of 
ACAD, where the x-axis is chosen to be perpendicular to the 
discontinuity, the surface curvature is such that T X >1, T y ~ 1 at 
the discontinuities inside the pupil and T X >1, T y ~l elsewhere. 
This yields damped chromatic oscillations at the remapped 
discontinuities and Fresnel oscillations at the edges of the 
pupil (see right panel of FIG. 10 ). Note that FIG. 10 was 
generated using a one dimensional toy model that assumes 
Eq. 23 is separable, e.g r =1, as described in [Pueyo et~al. 
(20 1 1 )] . In practice at the saddle points of the optical surfaces, 
near the junction of two spiders for instance, y x >l, ©^<-1 and 
thus our separable model does not guarantee than in the true 
2D case chromatic edge oscillations might not be locally 
amplified. However even near the saddle points ACAD pro- 
vides a strongly converging remapping in the direction per- 
pendicular to the discontinuity and very little diverging 
remapping in the other direction. As a consequence r x »l 
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and r is smaller than but close to one. We thus predict that 
even at these locations chromatic ringing will not be ampli- 
fied. Even though ACAD is based on pupil remapping, its 
diffraction properties are qualitatively very different from 
PIAA coronagraphs since edge ringing is not amplified 5 
beyond the Fresnel regime at the pupil edges, and is attenu- 
ated near the discontinuities. We conclude that in most cases 
ACAD operates in a regime where edge chromatic oscilla- 
tions are not larger than classical Fresnel oscillations, and 
sometimes actually smaller. As a consequence the chroma- to 
ticity of this ringing can be mitigated using standard tech- 
niques developed in the Fresnel regime, and we do not expect 
this phenomenon to be a major obstacle to ACAD broadband 
operations. 

5 .3 Diffraction Artifacts in ACAD are no Worse than Fresnel 15 
Ringing 

Results obtained when applying our approach to a geom- 
etry similar to JWST. We used two 3 cm DMs of 64 actuators 
separated by 1 m. Their maximal surface deformation is 1 . 1 
pm, well within the stroke limit of current DM technologies. 20 
The residual light in the corrected PSF follows the secondary 
support structures and can potentially be further cancelled by 
controlling the DMs using an image plane based cost func- 
tion, see FIG. 23 

We have established that the diffractive chromatic oscilla- 25 
tions introduced by the fact that DM1 and DM2 are not 
located in conjugate planes is no worse than classical Fresnel 
ringing from the aperture edges and can be mitigated using 
well-know techniques which have been developed for this 
regime. While a quantitive tradeoff study of how to design a 30 
high contrast instrument which minimizes such oscillations is 
beyond the scope of this paper, we briefly remind their quali- 
tative essence to the reader: 

The edges of the discontinuities in the telescope aperture 
can be smoothed via pupil apodization before DM1. 35 
This solution is not particularly appealing as it requires 
the introduction of a transmissive, and thus dispersive, 
component in the optical train. 

The distance between the two DMs can be reduced. Indeed 
the DMs deformations presented herein, for 3 cm DMs 40 
separated by 1 m, are all <1 pm while current technolo- 
gies allow for deformations of several microns. As the 
edge ringing scales as Z/D' chromatic oscillations will 
be reduced by decreasing Z. Since the DM surfaces scale 
as D 2 /Z reducing Z will increase the DM deformations 45 
but have little impact on the feasibility of our solution as 
current DM technologies can reach 4 pm strokes. 

The coronagrahic apodizer can be placed in a plane that is 
conjugate to the DM1. This can be achieved by re- 
imaging DM2 through a system of oversized optics (the 50 
over-sizing factor increases steeply when the pupil 
diameter decreases). By definition there are no Fresnel 
edge oscillations in such a plane. Alternatively a coro- 
nagraph without a pupil apodization (amplitude orphase 
mask in the image plane) can be used, and in this con- 55 
figuration it is only sufficient for the optics to be over- 
sized. 

Note that these three solutions are not mutually exclusive 
and that only a full diffractive analysis, which uses robust 
numerical propagators that have been developed based on Eq. 60 
23, can quantitatively address the trade-offs discussed above. 
The development of such propagators is our next priority. In 
[Pueyo et~al. (2011)] we laid out the theoretical foundations 
of such a numerical tool in the case of circularly symmetric 
pupil remapping and this solution has been since then prac- 65 
tically implemented, as reported by [Krist et~al, (2010)]. 
Generalizing this method to a tractable propagator in the case 


of arbitrary remapping is a yet unsolved computational prob- 
lem. In the meantime we emphasize that while the spectral 
bandwidth of coronagraphs whose incoming amplitude has 
been corrected using ACAD will certainly be limited by edge 
diffraction effects, but these effects are no worse than Fresnel 
ringing and can thus be mitigated using optical designs which 
are now routinely used in high contrast instruments (see 
[Vainaud et-al. (2010)] for such discussions). For the remain- 
der of this paper we thus assume the diffractive artifacts have 
been adequately mitigated and we compute our results 
assuming a geometric propagation between DM1 and DM2. 

When using two sequential DMs as a pure amplitude actua- 
tor (e.g., assuming a flat wavefront after the second DM) the 
relationship between the surface deformations and the outgo- 
ing electrical field is given by the Huygens integral. We 
derived, under the second order approximation for the Huy- 
gens integrand, a simplified expression of this diffractive 
transfer function: 
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The computational cost associated with the diffractive 
transfer function in the above equation and also Eq. (59) 
makes it almost impossible to devise a global inversion 
method that captures both the remapping of large amplitude 
errors and the diffractive properties of the dual DMs unit. 
Procedure to Include Diffraction 

As a consequence in order to find DM shapes that effi- 
ciently compensate for the secondary support structures, we 
provide the procedure outlined below: 

1 . First invert the problem in the ray optics approximation 
in order to map out the pupil imprint of the spiders as 
much as possible, as discussed in this specification. 

2. Second compute in the SR-approximation the diffractive 
response associated with the derived DM shapes using 
Eq. (59). Assess the broadband properties of the dual 
DM unit in the presence of diffractive effects. 

3. Using a small perturbation superposed to the deforma- 
tions calculated in step 1, compute, around the equilib- 
rium point set by step 1 , the sensitivity matrix for both 
DMs seen through the chosen coronagraph. This matrix 
is calculated using the diffractive models in step 2, e.g. in 
the SR-approximations. We conduct an iterative closed 
loop. 

Results for the Astronomical Telescope Application 

As an illustration we have applied the entire methodology 
described above to a geometry analogous to the AFTA tele- 
scopes. We have in particular successfully implemented a 
numerical integrator capable of capturing the diffractive 
effects in Eq. (59) that are the sources of the bandwidth limit 
in ACAD. Our results are shown on FIG. 26 where we calcu- 
lated, assuming that the DMs have been controlled according 
to step 1, the field in a plane conjugate with the telescope 
pupil. When introducing a Focal Plane Mask, this plane 
would be the Lyot plane of the coronagraph and its wavefront 
topology drives the contrast floor and bandwidth when using 
linear DM control loops. From a ray optics perspective 
ACAD “fills up as much as possible” the area obscured by the 
secondary support structure FIG. 3, while in the diffractive 
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regime we can see that ACAD replaces the deep and sharp 
amplitude errors from supports with shallow and smooth 
ripples. The chromatic nature of these oscillations sets the 
ultimate bandwidth of ACAD. We see that ACAD only 
weakly impacts the bandwidth of the starlight suppression 5 
system, as shown in FIG. 27 where we have combined these 
diffractive models with a Radial Apodization Vector Vortex 
Coronagraph (RA-WC, Mawet et al., in prep., and obtained 
contrasts as deep at 3x10 at $2\lambda/D$ over the hill range 
of wavelengths from 500 mn to 700 nm, a bandwidth of over to 
30%. The combination of FIG. 23 (APLC or PIAAC) and 
FIG. 27 illustrates the second advantage of ACAD: it is fun- 
damentally “coronagraph agnostic”. We have conducted dif- 
fractive simulations using the RA-WC since it is one of the 
most promising technologies for obscured apertures 70% 15 
throughput in its high transmissivity configuration and theo- 
retical IWA as close as 1 .5 X/D. Smaller I WA solutions should 
they be available later on. 

6 Results 

6.1 Application to Future Observatories 20 

6.1.1 JWST 

We have illustrated each step of the calculation of the DM 
shapes in §. 4 using a geometiy similar of JWST. This con- 
figuration is somewhat a conservative illustration of an on- 
axis segmented telescope as it features thick secondary sup- 25 
ports and a “small” number of segments whose gaps diffract 
light in regions of the image plane close to the optical axis (the 
first diffraction order of a six hexagons structure is located at 
-3X/D). In order to assess the performances of ACAD on such 
an observatory architecture we chose to use a coronagraph 30 
designed around a slightly oversized secondary obscuration 
of diameter 0.25 D, with a focal plane mask of diameter 3A./D, 
an IWA of 5X/D and an OWA of 30X/D. The field at the 
entrance of the coronagraph after remapping by the DMs is 
shown on the top right panel of FIG. 11 . The DM surfaces, 35 
calculated assuming 64 actuators across the pupil (N=64 in 
the Fourier expansion) and DMs of diameter 3 cm separated 
byZ=l m, are shown on the middle panel ofFIG. 11 . They are 
well within the stroke limit of current DM technologies. The 
surfaces were calculated by solving the reverse problem over 40 
an even grid of 10 cutoff low-pass spatial frequencies ranging 
between 30 and 70 cycles per apertures for the tapering kernel 
co. The value yielding the best contrast was chosen. Note that 
the optimal cutoff frequency depends on the spatial scale of 
the discontinuities, and that higher contrasts could be 45 
obtained by choosing a set of two convolution kernels for the 
reverse problem and finding the optimal solution using a finer 
grid. However, the results in the bottom row of FIG. 11 are 
extremely promising. FIG. 12 shows a contrast improvement 
of a factor of 100 when compared to the raw coronagraphic 50 
PSF, which is quite remarkable for an algorithm which is not 
based on an image-plane metric. These results illustrate that 
even with a very unfriendly aperture similar to JWST one can 
obtain contrasts as high as envisioned for upcoming Ex -AO 
instruments, which have been designed for much friendlier 55 
apertures. While we certainly do not advocate to use such a 
technique on JWST, this demonstrates that ACAD is a pow- 
erful tool for coronagraphy with on-axis segmented aper- 
tures. 

6.1.2 Extremely Large Telescopes 60 

PSFs resulting from ACAD when varying the number and 

thickness of secondary support structures while maintaining 
their covered surface constant (FIG. 15 ). The surface area 
covered in this example is 5096 greater than in the TMT 
example shown on FIG. 13 . As the spiders get thinner their 65 
impact on raw contrast becomes smaller and the starlight 
suppression after DM correction becomes bigger. For a rela- 


36 

tively small number of spiders (<12) the contrast improve- 
ment on each single structure is the dominant phenomenon, 
regardless of the number of spiders. ELTs designed with a 
moderate to large number of thin secondary support structure 
(6 to 1 2) present aperture discontinuities which are easy to 
correct with ACAD. 

We now discuss the case of Extremely Large Telescopes 
and provide an illustration using the example of the Thirty 
Meter Telescope. We considered the aperture geometry 
shown on the top left panel ofFIG. 13 . It consists of a pupil 37 
segments across in the longest direction and a secondary of 
diameter ~0. 1 2D which is held by three main thick struts and 
six thin cables. As seen on the bottom left panel ofFIG. 13 the 
impact of segment gaps is minor as they diffract light beyond 
the OWA of the coronagraph. When using a coronagraph with 
a laiger OWA the segment gaps will have to be taken into 
account, and will have to be mitigated using DMs with a 
larger number of actuators. In order to obtain first order esti- 
mates of the performances of ACAD on the aperture geom- 
etry shown on the top left panel ofFIG. 13 , we chose to use a 
coronagraph designed around a slightly oversized secondary 
obscuration of diameter 0. 1 5 D, with a focal plane mask of 
6A7D diameter, an IWA of 4X/D and an OWA of 30A/D. The 
field at the entrance of the coronagraph after remapping by the 
DMs is shown on the top right panel of FIG. 13 . The DM 
surfaces, calculated assuming 64 actuators across the pupil 
(N=64 in the Fourier expansion) and DMs of diameter 3 cm 
separated by Z=1 m, are shown on the middle panel ofFIG. 
13 . They are well within the stroke limit of current DM 
technologies. The surfaces were calculated by solving the 
reverse problem over an even grid of 10 cutoff low -pass 
spatial frequencies ranging between 30 and 70 cycles per 
apertures for the tapering kernel w. The value yielding the 
best contrast was chosen. The final PSF is shown on the 
bottom right panel ofFIG. 13 and features a high contrast dark 
hole with residual diffracted light at the location of the spi- 
ders’ diffraction structures. The impact on coronagraphic 
contrast of secondary supports was thoroughly studied by 
[Martinez et~al. (2008)]. They concluded that under a 90% 
Strehl ratio, the contrast in most types of coronagraphs is 
driven by the secondary support structures to levels ranging 
from 10 -4 to 10 -5 . This, in turn, leads to a final contrast after 
post-processing (called Differential Imaging) of ~10 _7 -10 -8 . 
FIG. 14 shows that using ACAD on an ELT pupil yields 
contrasts before any post -processing which are comparable to 
the ones obtained by [Martinez et~al. (2008)] after Differen- 
tial Imaging. This demonstrates that should two sequential 
DMs be integrated into a future planet finding instrument, 
setting their surface deformation according to the methodol- 
ogy presenting above would allow this instrument to perform 
its scientific program at a very high contrast. Moreover the 
surface of the DMs could be adjusted to mitigate for the effect 
of missing segments at the surface of the primary (when for 
instance the telescope is operating while some segments are 
being serviced). 

6.2 Hypothetical Cases 

6.2.1 Constant Area Covered by the Secondary Support 
Structures 

In the case of ELTs with large number of small segments 
(when compared to the aperture size), gaps diffract light far 
from the optical axis (see FIG. 13 for an example). The 
secondary support structures are then the major source of 
unfriendly coronagraphic diffracted light. Under the assump- 
tion that thick structures are necessary to support the heavy 
secondary over the very large ELT pupils, one can use the 
aperture area covered by the spiders as a proxy of the second- 
ary lift constraint. We have thus explored a series of geom- 
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etries for which the number of spiders increases as they get 
thinner while the overall area covered by the secondary sup- 
port structures remains constant. In the examples shown from 
FIG. 15 to FIG. 18 , the area covered by the secondary support 
structures is 1.5 times greater than in the TMT geometry 5 
discussed above. In all cases we used a coronagraph with a 
central obscuration of 0. 1 5 D, with a focal plane mask of 6>JD 
diameter, and IWA of 4X/D and an OWA of 30X/D. The 
surfaces were calculated by solving the reverse problem over 
an even grid of 10 cutoff low-pass spatial frequencies ranging to 
between 30 and 70 cycles per apertures for the tapering ker- 
nel. The value yielding the best contrast was chosen. This 
exercise leads to several conclusions pertaining to the perfor- 
mances of ACAD with various potential ELT geometries. 

Clocking of the Spiders with Respect to the DM 15 

The top two panels of FIG. 15 illustrate the importance of 
the clocking of the spiders with respect to the DMs actuator 
grid (or the Fourier grid in our case). When the secondary 
support structures are clocked by 45° with respect to the DM 
actuators they are much more attenuated by ACAD, thus 20 
yielding higher contrast. This is an artifact of the Fourier basis 
set chosen and would be mitigated by using DMs whose 
actuator placement presents circular and azimuthal symme- 
tries [Watanabe et~al. (2008)]. 

Annulus in the PSF with a Large Number of Spiders 25 

When the number of secondary support struts becomes 
very large (>20), an interesting phenomena occurs in the raw 
PSF: the spiders diffract light outside an annulus of radius 
Kspirferj/Jt^/D, just as spiderweb masks do in the case of 
shaped pupil coronagraphs [Vanderbei et~al. (2003)]. The 30 
“bump” located beyond that spatial frequency is more diffi- 
cult to attenuate using the DMs (see FIG. 15 for an illustra- 
tion). ACAD creates small ripples at the edges of the 
remapped discontinuities and when too many discontinuities 
are in the vicinity of each other, then these ripples interfere 35 
constructively and hamper the starlight extinction level 
yielded by ACAD. 

A Lot of Thin Spiders is More Favorable Than a Few Thick 
Spiders 

In general decreasing the width of the spiders while 40 
increasing their number is beneficial to the contrast obtained 
after ACAD as illustrated on the radial averages on FIG. 16 
and FIG. 18 . When one increases the number of spiders while 
decreasing their width in a classical coronagraph, the peak 
intensity of the diffraction pattern of one spider decreases as 45 
the squared width of the spider. The radially averaged contrast 
improvement without ACAD is then somewhat lesser than the 
square of the spider thinning factor as it is mitigated by the 
increasing number of spiders. When using ACAD the spiders 
are seen by the coronagraph as much thimier than they actu- 50 
ally are (by a factor x) and thus the peak intensity of their 
diffraction pattern is lower by a factor of x 2 , Our numerical 
experiments show that x increases when the spider width 
decreases. As consequence, the overall contrast gain after 
ACAD when decreasing the width of the spiders while 55 
increasing their number is greater than in the case of a clas- 
sical coronagraph. When designing ELT secondary support 
structures and planning to correct for them using ACAD, 
increasing the number of spiders to 8 or even 12 has a ben- 
eficial impact on contrast as it enables each discontinuity to 60 
become thinner and thus to be corrected to higher contrast 
using the DMs. The PSFs of apertures with more than 12 
spiders present diffraction structures which are poorly suited 
for correction with square DMs. While the contrast resulting 
from applying ACAD to such apertures is still a decreasing 65 
function of the number of spiders, FIG. 18 shows that the net 
contrast gain brought by the DM based remapping is smaller 
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than in the more gentle cases of 12 spiders. The study pre- 
sented on FIG. 15 to FIG. 18 remains to be fully optimized for 
each potential design of an ELT planet finding instrument (in 
particular using a finer grid of cutoff spatial frequencies). It 
however demonstrates the flexibility of ACAD for various 
aperture geometries and provides a first order rule of thumb to 
design telescope apertures which are friendly to direct imag- 
ing of exo-planets: “A lot of thin spiders is more favorable 
than a few thick spiders”. In practice the number of spiders 
will be limited by effects not treated in our analysis such as the 
mechanical rigidity, requirements on the perfection of their 
periodic spacing and glancing reflections from the sides of 
multiple spiders. We thus advocate that, should future ELTs 
be built with high contrast exo-planetary science as a main 
scientific driver, then such effect ought to be thoroughly ana- 
lyzed as a large numbers of thin spiders is more favorable 
from a contrast standpoint when using ACAD to mitigate for 
pupil amplitude asymmetries. 

6.2.2 Monolithic on- Axis Apertures. 

PSFs resulting from ACAD when varying the number and 
thickness of secondary support structures. As the spiders get 
thimier their impact on raw contrast becomes lesser and the 
starlight suppression after DM correction becomes greater. In 
this case ui was optimized on a very fine grid and the aperture 
we clocked in a favorable direction with respect to the Fourier 
basis. 

When discussing the case of JWST we stressed the com- 
plexity associated with the optimization of ACAD in the 
presence of aperture discontinuities of varying width. Carry- 
ing out such an exercise would be extremely valuable to study 
the feasibility of the direct imaging of exo-earths with an 
on-axis segmented future flagship observatory such as 
ATLAST [Postman et~al. (2010)]. However, such an effort is 
computationally heavy and thus beyond the scope of the 
present paper, which focuses on introducing the ACAD meth- 
odology and illustrating using key basic examples. 

So far, none of the examples in this manuscript demon- 
strate that ACAD can yield corrections all the way down to the 
theoretical contrast floor that is set by the coronagraph design. 
When seeking to image exo-earths from space, future mis- 
sions will need to reach this limit. In order to explore this 
regime, we conducted a detailed study of a hypothetical on- 
axis monolithic telescope with four secondary support struts. 
To establish the true contrast limits we varied the thickness of 
the spiders for each geometry. The surfaces were calculated 
by solving the reverse problem over an even grid of 70 cutoff 
low-pass spatial frequencies ranging between 30 and 70 
cycles per apertures for the tapering kernel. The value yield- 
ing the best contrast was chosen. In all cases we used a 
coronagraph with a central obscuration of 0. 1 5 D, with a focal 
plane mask of 677D diameter, and IWA of 4X/D, and an OWA 
of 30>JD. Note that when using coronagraphs relying on 
pupil apodization these results can be readily generalized to 
larger circular secondary obscurations, at a loss in IWA (as 
shown on FIG. 2 ). Moreover we clocked the telescope aper- 
ture by 45° with respect to the grid of Fourier modes. We 
found that, indeed, the theoretical contrast floor set by the 
coronagraph design is met for thin spiders (0.02 D) and it is 
very close to being met for spiders only half the thickness of 
the ones currently equipping the Hubble Space Telecope 
(0.05 D), see FIG. 19 and FIG. 20 . Even in the case of thick 
stmts (0.1 D) we find contrasts an order of magnitude higher 
than in the similar configuration on the top panel of FIG. 10, 
due to our thorough optimization of the cutoff frequency of 
the tapering kernel and careful clocking of the aperture with 
respect to the actuators. On-axis telescopes are thus a viable 
option to image earth-analogs from space: their secondary 
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support structures can be corrected down to contrast levels 
comparable to the target contrast of recent missions concept 
studies [Guyon et~al. (2008), Trauger et~al. (2010)]. Since 
the baseline wavefront control architecture for future space 
coronagraphs relies on two sequential DMs, ACAD does not 5 
add any extra complexity to such missions and merely con- 
sists of controlling the DMs in order to optimally compensate 
for the effects of asymmetric aperture discontinuities. 

7 Discussion 

7.1 Field Dependent Distortion 10 

Because ACAD relies on deforming the DMs’ surfaces in 

an aspherical fashion, off-axis wavefronts seen through the 
two DMs apparatus will be distorted, just as in a PI AA coro- 
nagraph [Martinache et~al. (2006)]. However the asphericity 
of the surfaces in the case of ACAD operating on reasonably 15 
thin discontinuities, is much smaller than in a PIAA remap- 
ping unit. FIG. 21 shows the impact on off-axis PSFs of such 
a distortion in the worse-case scenario of a geometry similar 
to JWST. We demonstrate that most of the flux remains in the 
central disk of radius X/D for all sources in the field of view of 20 
the coronagraphs considered here (all the way to 30 cycles per 
capture). We conclude that, because of the small deforma- 
tions of the DMs, PSF distortion will not be a major hindrance 
in exo-planet imaging instruments whose DMs are controlled 
in order to mitigate for discontinuities in the aperture. 25 

7.2 Impact of Wavefront Discontinuities in Segmented Tele- 
scopes. 

7.2.1 General Equations in the Presence of Incident Phase 
Errors and Discontinuities 

So far we have treated primary mirrors’ segmentation as a 30 
pure amplitude effect. In reality the contrast floor in seg- 
mented telescopes will be driven by both phase and amplitude 
discontinuities; here we explore the impact of phase errors 


amplitude at DM2. In that case treatments of these phenom- 
enons that have already been discussed in the literature for 
monolithic apertures are still valid for small enough phase 
errors and smooth enough remapping functions. For ACAD 
remapping this smoothness constraint is naturally enforced 
by the limited number of actuators across the DM surface. In 
this case phase to amplitude conversioncan in principle be 
corrected using DM1. (2) Remapped phase discontinuities 
can be corrected for a finite number of spatial frequencies 
using a continuous phase sheet deformable mirror. We illus- 
trate this partial correction over a 20% bandwidth using 
numerical simulations of a post- ACAD half dark hole created 
by superposing a small perturbation, computed using a linear 
wavefront control algorithm, to the ACAD DM2 surface. 

If the incoming wavefront is written as Alq and the solution 
of the Monge- Ampere Equation for DM1 as lq 0 then one can 
conduct the analysis in Eq. 9 to 19 using h^h^+Alq. Under 
the assumption that surface of DM2 is set as h 2 in order to 
conserve Optical Path Length then one can re-write the 
remapping as (fj.jq) defined by; 


a/j, 


dx 


Vi(x2’y2)’gi( x 2>y2)) 


fi(x 2 ,y 2 )-x 2 

z 


( 63 ) 


a/j, 


dy 


V 1 ( r 2’ y 2 '■gl^2->'2 )l 


j i (■*; ■ y; ) ~ yi 

Z 


Moreover if edge ringing has been properly mitigated then 
the ray optics solution is valid and the field at DM2 can be 
written as: 
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and discontinuities occurring before two DMs whose sur- 
faces have been set using ACAD. There are two main phe- 
nomena to be considered. The first is the conversion of the 
incident wavefront phase before 


DM1: —A/i, 

A 

into amplitude at the second mirror. The second is the projec- 
tion of this wavefront phase into a remapped phase errors at 


DM 2: — A/ii(/i(.V!. y 2 ), gi(xi,y 2 )). 

Since the remapping unit is designed using deformable mir- 
rors, both DM1 and DM2 a complete correction could be 60 
attained in principle. However, the deformable mirrors are 
continuous while Alq presented discontinuities. Thus, com- 
plete corrections for segmented mirrors might not be 
achieved in practice. Below we discuss the following two 
main points. (1) Even if the phase wavefront error Alq has 65 
discontinuities, the phase errors within in segment still drive 
the phase to amplitude conversion and thus the propagated 


7.2.2 Impact on the Amplitude after ACAD 

Broadband wavefront correction (20% bandwidth around 
700 nm) with a single DM in segmented telescope with dis- 
continuous surface errors. Top Left: wavefront before the 
coronagraph. Top Right: broadband aberrated PSF with DM 
at rest. Bottom Left: DM surface resulting from the wavefront 
control algorithm. Bottom Right: broadband corrected PSF. 
Note that the wavefront control algorithm seeks to compen- 
sate for the diffractive artifacts associated with the secondary 
support structures: it attenuates them on the right side of the 
PSF while it strengthens them on the left side of the PSF. As 
a result the DM surface becomes too large at the pupil spider’ s 
location and the quasi-linear wavefront control algorithm 
eventually diverges. 

We first consider the amplitude profile in Eq. 64. it is 
composed of two factors: the remapped telescope aperture, 
E iwl (I’i,gi), and the determinant of Id+7f [hj. 

The first condition necessary for the incoming wavefront 
not to perturb the ACAD solution is: Ah, is such that the 
remapping is not modified at the pupil locations where the 
telescope aperture is not zero E flM /0. This results in the 
following conditions: 
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At the locations where E^^^ there is no light illuminat- 
ing the discontinuous wavefront and thus the large local 
slopes at these location have no impact on the remapping 10 
functions (f 1: gj. These conditions are not true in segmented 
telescopes that are not properly phased, for which the tip-tilt 
error over each segment can reach several waves. However 
under the assumption that the primary has been properly 
phased (for instance the residual rms wavefront after phasing 1 5 
is expected to be ~Vio th of a wave, similar to values expected 
for JWST NIRCAM) then these conditions are true within the 
boundaries of each segment. Moreover, while the local wave- 
front slopes at the segment’s discontinuities do not respect 
this condition the incident amplitude at these points is ^i DMl 20 
(X, y)=0 and they thus do not perturb the ACAD remapping 
solution. 

The second necessary condition resides in the fact that the 
determinant of Id+ H (hj is not equal to det[Id+ (h, 0 )] at 
the pupil locations where the telescope aperture is not zero 25 
ought not have a severe impact on contrast. One can use the 
linearization in [Loeper Rapetti (2005)] to show that: 
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The perturbation term AA(Ah t ) corresponds to the full 
non-linear expression of the phase to amplitude conversion of 
wavefront errors that occurs in pupil remapping units. In 
[Pueyo et~al. (2011)] we derived a similar expression in the 45 
linear case, when Ah 1 «X and showed that in the pupil 
regions where the beam is converging this phase to amplitude 
conversion was enhanced when compared to the case of a 
Fresnel propagation. In a recent study [Krist et~al. (2011)] 
presented simulations predicting that this effect was quite 50 
severe in PIAA coronagraphs and can limit the broadband 
contrast after wavefront control unless DMs were placed 
before the remapping unit. In principle ACAD will not suffer 
from this limitation as the first aspherical surface of the 
remapping unit is actually a Defonnable Mirror that can 55 
actually compensate for All,, before any phase to amplitude 
wavefront modulation occurs. Devising a wavefront control- 
ler that relies on DM1 requires moreover a computationally 
efficient model to propagate arbitrary wavefronts through 
ACAD. Such a tool was developed in [Krist et~al. (2010)] 60 
assuming azimuthally symmetric geometries. Since devising 
such a tool in ACAD’ s case, in the asymmetric case, repre- 
sents a substantial effort well beyond the problem of prescrib- 
ing ACAD DM shapes, we chose not to include such simula- 
tions in the present manuscript. Since we are using 65 
Deformable Mirrors with a limited number of actuators, 
ACAD remapping is in general less severe than in the case of 
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PIAA. We thus expect the results regarding the wavefront 
correction before the remapping unit reported in [Krist et~al. 
(2011)] to hold. This is provided that the DM actuators can 
adequately capture the high spatial frequency content of Afq 
to create a dark hole in the coronagraphic PSF. We tackle this 
particular aspect next when discussing the case of phase 
errors, in the absence of wavefront phase to amplitude con- 
version. Once again note that while the local wavefront cur- 
vatures are very large at the segment’s discontinuities, the 
incident amplitude at these points is E DA/1 (x,y)=0 and thus 
they do not have an impact on the ACAD phase to amplitude 
modulation. In practice if the DM is not exactly located at a 
location conjugate to the telescope pupil the actual wavefront 
discontinuities will be slightly illuminated and mi gilt perturb 
the remapping functions and the phase to amplitude conver- 
sion. While this might tighten requirements regarding the 
positioning of DM1 in the direction of the optical axis we do 
not expect this effect to be a major obstacle to successful 
ACAD implementations. 

7.2.3 Impact on the Phase after ACAD 

In practice, when Alij presents discontinuities, the surface 
of DM2 cannot be set to the deformation h 2 that conserves 
Optical Path Length, since we work under the assumption that 
the DM has a continuous phase-sheet. While this has no 
impact on the discussion above regarding the amplitude of 
E/W 2 - since DM1 is solely responsible for this part, it ought to 
be taken into account when discussing the phase at DM2. 
Under the assumption that Alq does not perturb the nominal 
ACAD remapping function then one can show that the phase 
at DM2 is: 


arg[£ D M:U-;, >2) 1 == y (A/ifi/iW >2), <?i(*2. A'z)) + A/i 2 (jc 2 , >2)) 

where Ah 2 (x 2 , y 2 ) is a small continuous surface deformation 
superposed to the ACAD shape of DM2 and Ah 1 (f 1 °(x 2 , y 2 ), 
gi°(x 2 , y 2 )) is the telescope OPD seen through the DM based 
remapping unit. This second term presents phase discontinui- 
ties whose spatial scale has been contracted by ACAD. When 
these discontinuities are very small then their high spatial 
frequency content does not disrupt the ability of DM2 to 
correct for low to mid-spatial frequency wavefront errors 
wavefront errors. However as the discontinuities become 
larger their high spatial frequency content can fold into the 
region of the PSF that the DMs seek to cancel. These “fre- 
quency folding” speckles are highly chromatic [Give’ on 
et~al. (2006)] and can have a severe impact on the spectral 
bandwidth of a coronagraph whose wavefront is corrected 
using a continuous DM. 

In order to assess the impact of this phenomenon, we 
conducted a series of simulations based on single DM wave- 
front control algorithm that seeks to create a dark hole in one 
half of the image plane at in as in [Borda&Traub (2006)]. We 
use the example of a geometry similar to JWST and work 
under the assumption that the discontinuous wavefront inci- 
dent to the coronagraph has the same spatial frequency con- 
tent as a JWST NIRCAM Optical Path Difference that has 
been adjusted to 70 nm mis in order to mimic a visible Strehl 
similar to the near-infrared Strehl of JWST. The non-linear 
wavefront and sensing and control problem associated with 
phasing a primary mirror to such level of precision is 
undoubtedly a colossal endeavor and is well beyond the scope 
of this paper. In this section we work under the assumption 
that the primary mirror either has been phased to such a level, 
that the wavefront discontinuities are no larger than 200 mn 
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peak to valley or that the wavefront has been otherwise cor- 
rected down to this specification using a segmented Deform- 
able Mirror that is conjugate with the primary mirror. More- 
over we assume (1) that the residual post-phasing wavefront 
map has been characterized and can be used in order to build 
the linear model underlying the wavefront controller (2) the 
focal plane wavefront estimator (carried using DM diversity 
as in [Borda. & Traub (2006)] for instance) is capable to yield 
an exact estimate of the complex electrical field at the science 
camera. Underlying this last assumption is the overly opti- 
mistic premise that wavefront will remain unchanged over the 
course of each high-contrast exposure. While this is not a 
realistic assumption one could envision the introduction of 
specific wavefront sensing schemes, with architectures simi- 
lar to the one currently considered for low order wavefront 
sensors on monolithic apertures [Guyon et~al. (2009), Wal- 
lace et~al. (2011)], or using a separate metrology system. The 
results presented here are thus limited to configurations for 
which segment phasing will be dynamically compensated 
using specific sensing and control beyond the scope of this 
paper. As this section merely seeks to address the controlla- 
bility of wavefront errors in segmented telescopes we chose to 
conduct our simulations with a perfect estimator. Finally we 
use the stroke minimization wavefront control algorithm pre- 
sented in [Pueyo et~al. (2009)] to ensure convergence for as 
many iteration as possible. We first tested the case of a seg- 
mented telescope in the absence of ACAD, using a azimuth- 
ally symmetric coronagraph and a single DM. We sough to 
create a Dark Hole between 5 and 28 X 0 /D under a 20% 
bandwidth with X o =700 lun. FIG. 22 shows the results of such 
a simulation. The DM can indeed correct for the discontinui- 
ties over a broadband in one half of the image plane. However 
the wavefront control algorithm seeks to compensate for the 
diffractive artifacts associated with the secondary support 
structures: it attenuates them on the right side of the PSF 
while it strengthens them on the left side of the PSF. As a 
result the DM surface becomes too large at the pupil spider’s 
location and the quasi-linear wavefront control algorithm 
eventually diverges for contrasts ~10 6 . 

Broadband wavefront correction (20% bandwidth around 
700 nm) in a segmented telescope whose pupil has been 
re-arranged using ACAD. The surface of the first DM is set 
according to the ACAD equations. The surface of the second 
DM is the sum of the ACAD solution and a small perturbation 
calculated using a quasi-linear wavefront control algorithm. 
Top Left: wavefront before the coronagraph. Note that the 
ACAD remapping has compressed the wavefront errors near 
the struts and the segment gaps. Top Right: broadband aber- 
rated PSF with DMs set to the ACAD solution. Bottom Left: 
perturbation of DM2’ s surface resulting from the wavefront 
control algorithm. Bottom Right: broadband corrected PSF. 
The wavefront control algorithm now yields a DM surface 
that does not feature prominent deformations at the location 
of the spiders. Most of the DM stroke is located at the edge of 
the segments, at location of the wavefront discontinuities. 
There, the DM surface eventually becomes too large and the 
quasi-linear wavefront control algorithm diverges. However 
this occurs higher contrasts than in the absence of ACAD. 

We then proceeded to simulate the same configuration in 
the presence of two DMs whose surface at rest was calculated 
using ACAD. Since there does not exist a model yet to propa- 
gate arbitrary wavefronts through ACAD (the models in 
[Krist et~al. (2011)] only operate under the assumption of an 
azimuthally symmetric remapping) we can only use the sec- 
ond DM for wavefront control. We work under the assump- 
tion that the incident wavefront does not perturb the nominal 
ACAD remapping (which is true in the case of the surface 
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map we chose for our example) and that the arguments in 
[Krist et~al. (2010)] hold so that phase to amplitude conver- 
sion in ACAD can be compensated by actuating DM1 . Fre- 
quency folding will then be the phenomenon responsible for 
5 the true contrast limit. In this section we are interested in 
exploring how this impacts the controllability of wavefront 
discontinuities using continuous phase-sheet DMs. We used a 
azimuthally symmetric coronagraph and superposed our 
wavefront control solution to DM2. We sought to create a 
to Dark Hole between 5 and 23 >^ 0 /D under a 20% bandwidth 
with Xq= 700 nm. FIG. 22 shows the results of such a simula- 
tion. When the incident wavefront is small enough it is indeed 
possible to superpose a “classical linear wavefront control” 
solution to the non-linear ACAD DM shapes in order to carve 
15 PSF dark holes. The wavefront control algorithm now yields 
a DM surface that does not feature prominent deformations at 
the location of the spiders. Most of the DM stroke is located 
at the edge of the segments, at location of the wavefront 
discontinuities and seek to correct the frequency folding 
20 terms associated with such discontinuities . At these locations 
the DM surface eventually becomes too large and the linear 
wavefront control algorithm diverges. However this diver- 
gence occurs at contrast levels much higher than when the 
ACAD solution is not applied to the DMs. These simulations 
25 show that indeed discontinuous phases can be corrected using 
the second DM of a ACAD whose surfaces have preliminary 
been set to mitigate the effects of spiders and segment gaps. 
7.3 Ultimate Contrast Limits 

Assuming that edge ringing has been properly mitigated so 
30 that the ray optics approximation underlying the calculation 
of the DMs shapes is valid, one can wonder about the ultimate 
contrast limitations of the results presented in this manu- 
script. Increasing the number of actuators would have dra- 
matic effects on contrast if the actuator count would be such 
35 that N>D/d where d is the scale of the aperture discontinui- 
ties. Unfortunately current DM technologies are currently far 
from such a requirement and the solutions presented here are 
in the regime where N«D/d. In this regime N only has a 
marginal influence on contrast when compared to the impact 
40 of the cutoff frequency of the tapering kernel. In the regime 
described here varying the actuator count only changes the 
size of the corrected region. 

The residual PSF artifacts in FIGS. 11 to 20 follow the 
direction of the initial diffraction pattern associated with sec- 
45 ondary support structures and segments. When addressing the 
problem of aperture discontinuities by solving the Monge- 
Ampere Equation, ACAD calculates the DM shapes based on 
a pupil plane metric and thus mostly focuses on attenuating 
these stnictures with little regard to the final contrast. It is 
50 actually quite remarkable that such a pupil-only approach 
yields levels of starlight extinction of two to three orders of 
magnitude. A more appropriate metric would be the final 
intensity distribution in the post-coronagraphic image plane. 
However, as discussed in §. 3 classical wavefront control 
55 algorithms based on a linearization of the DMs deformations 
around local equilibrium shapes (such as the ones presented 
in [Bora & Traub (2006), Give’on et~al. (2007)] in the one 
DMcaseor [Pueyo et~al. (2009)] foroneortwo DMs)camiot 
be used to compensate the full aperture discontinuities. This 
60 is illustrated in FIG. 22 where the DM surface in the vicinity 
of spiders becomes too laige after a certain number of itera- 
tions, which leads the iterative algorithm to diverge. When 
attempting to circumvent this problem by recomputing the 
linearization at each iteration, we managed to somewhat sta- 
65 bilize the problem for a few iterations and reached marginal 
contrast improvements, but the overall algorithm remained 
unstable unless a prohibitively small step size was used. This 
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is the problem which motivated our effort to calculate the DM 
shapes as the full non-linear solution of the Monge- Ampere 
Equation. While doing so yields significant contrast improve- 
ments in both the case of JWST like geometries, TMT and on 
axis-monolithic apertures similar, this approach does not give 
a proper weight to the spatial frequencies of interest for high 
contrast imaging. We mitigated this effect by giving a strong 
weight to the spatial frequencies of interest (in the Dark Hole) 
when solving the Monge- Ampere Equation. 

The next natural step is thus to use non-linear solutions 
presented herein to correct for the bulk of the aperture dis- 
continuities and to serve as a starting point for classical lin- 
earized waveform control algorithms, as illustrated on FIG. 
24 . FIG. 25 indeed illustrates that when superposing an image 
plane based wavefront controller to the Monge-Ampere 
ACAD solution, the contrast can be improved beyond the 
floor shown on FIG. 12 . However, one DM solutions are of 
limited interest as they only operate efficiently over a finite 
bandwidth and over half of the image plane. ACAD yields a 
true broadband solution, and consequently it would be pref- 
erable to use the two DMs in the quasi-linear regime to 
quantify the true contrast limits of ACAD. In such a scheme 
the DM surfaces are first evaluated as the solution of the 
Monge-Ampere Equation and then adjusted using the image 
plane based wavefront control algorithm presented in [Pueyo 
et~al. (2009)]. However such an exercise requires efficient 
and robust numerical algorithms to evaluate Eq. 23 . Such tool 
only exist so far in the case of azimuthally symmetric remap- 
ping units [Krist et~al. (2010)]. Developing such numerical 
tools is thus of primary interest to both quantifying the chro- 
maticity and the true contrast limits achievable with on-axis 
and/or segmented telescopes. 

8 Conclusion 

We have introduced a technique that takes advantage of the 
presence of Deformable Mirrors in modem high-contrast 
coronagraph to compensate for amplitude discontinuities in 
on-axis and/or segmented telescopes. Our calculations pre- 
dict that this high throughput class of solutions operates under 
broadband illumination even in the presence of reasonably 
small wavefront errors and discontinuities. Our approach 
relies on controlling two sequential Deformable Mirrors in a 
non-linear regime yet unexplored in the field of high-contrast 
imaging. Indeed the mirrors’ shapes are calculated as the 
solution of the two-dimensional pupil remapping problem, 
which can be expressed as a non-linear partial differential 
equation called Monge-Ampere Equation. We called this 
technique Active Compensation of Aperture Discontinuities. 
While we illustrated the efficiency of ACAD using Apodized 
Pupil Lyot and Phase Induced Amplitude Coronagraph, it is 
applicable to all types of coronagraphs and thus enables one 
to translate the past decade of investigation in coronagraphy 
with unobscured monolithic apertures to a much wider class 
of telescope architectures. Because ACAD consists of a 
simple remapping of the telescope pupil, it is a true broadband 
solution. Provided that the coronagraph chosen operates 
under a broadband illumination, ACAD allows high contrast 
observations over a large spectral bandwidth as pupil remap- 
ping is an achromatic phenomenon. We showed that wave- 
length edge diffraction artifacts, which are the source of spec- 
tral bandwidth limits in PIAA coronagraphs (also based on 
pupil remapping), are no largerthan classical Fresnel ringing. 
We thus argued that they will only marginally impact the 
spectral bandwidth of a coronagraph whose input beam has 
been corrected with ACAD. The mirror deformations we find 
can be achieved, both in curvature and in stroke, with tech- 
nologies currently used in Ex- AO ground based instruments 
and in various testbeds aimed at demonstrating high-contrast 
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for space based applications. Implementing ACAD on a given 
on-axis and/or segmented thus does not require substantial 
technology development of critical components. 

For geometries analogous to JWST we have demonstrated 
5 that ACAD can achieve at least contrast ~10 -7 , provided that 
dynamic high precision segment phasing canbe achieved. For 
TMT and ELT, ACAD can achieve at least contrasts ~10 -s . 
For on-axis monolithic observatories the design contrast of 
the coronagraph can be reached with ACAD when the sec- 
10 ondary support structures are 5 times thinner than on HST. 
When they arejust as thick as HST contrasts as high as 10 s can 
be reached. These numbers are, however, conservative: an 
optimal solution can be obtained by fine tuning the control 
15 tenn in the Monge-Ampere Equation to the characteristic 
scale of each discontinuity. As our goal was to introduce this 
technique to the astronomical community and emphasize its 
broad appeal to a wide class of architectures (JWST,AT- 
LAST,HST,TMT,E-ELT) we left this observatory specific 
20 exercise for future applications. 

The true contrast limitation of ACAD resides in the fact that 
the Deformable Mirrors are controlled using a pre-corona- 
graph pupil based metric. However, as illustrated in FIG. 16 , 
the solution provided by ACAD can be used as the starting 
25 point for classical linearized waveform control algorithms 
based in image plane diagnostics. In such a control strategy, 
the surfaces are first evaluated as the solution of the Monge- 
Ampere Equation and then adjusted using the quasi-linear 
method presented in [Pueyo et~al . (2009)] . This control strat- 
30 egy requires efficient and robust numerical algorithms to 
evaluate the hill diffractive propagation in the remapped 
Fresnel regime. All the contrasts reported here are achieved 
without aberrations, and we showed that in practice, quasi- 
linear DM controls based on images at the science camera 
35 will have to be superposed to the ACAD solutions. Finally, as 
ACAD is broadly applicable to all types of coronagraphs, the 
remapped pupil can be used as the entry point to relax the 
design of coronagraphs that do operate on segmented aper- 
tures such as discussed in [Carlotti et~al. (2011), Guyon 
40 et~al. (2010)], also illustrated in FIG. 16 . ACAD is thus a 
promising tool for future high contrast imaging instruments 
on a wide range of observatories as it will allow astronomers 
to devise high throughput broadband solutions for a variety of 
coronagraphs. It only relies on hardware (Deformable Mir- 
45 rors) that have been extensively tested over the past ten years . 
Finally since ACAD can operate with all type of corona- 
graphs and it renders the last decade of research on high- 
contrast imaging technologies with off-axis unobscured aper- 
tures applicable to much broader range of telescope 
50 architectures. 
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The embodiments illustrated and discussed in this specifi- 
cation are intended only to teach those skilled in the art how 
to make and use the invention. In describing embodiments of 
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the invention, specific terminology is employed for the sake 
of clarity. However, the invention is not intended to be limited 
to the specific terminology so selected. The above-described 
embodiments of the invention may be modified or varied, 
5 without departing from the invention, as appreciated by those 
skilled in the art in light of the above teachings. It is therefore 
to be understood that, within the scope of the claims and their 
equivalents, the invention may be practiced otherwise titan as 

specifically described, 
to 

We claim: 

1 . An active optical beam shaping system, comprising: 
a first deformable mirror arranged to at least partially inter- 
15 cept an entrance beam of light and to provide a first 
reflected beam of light; 

a second deformable mirror arranged to at least partially 
intercept said first reflected beam of light from said first 
deformable mirror and to provide a second reflected 
20 beam of light; and 

a signal processing and control system configured to com- 
municate with said first and second deformable mirrors, 
wherein said signal processing and control system is con- 
figured to provide control signals to said first deformable 
25 mirror so as to configure a reflecting surface of said first 
deformable mirror to substantially conform to a calcu- 
lated surface shape, 

wherein said signal processing and control system is fur- 
ther configured to provide control signals to said second 
30 deformable mirror to configure a reflecting surface of 
said second deformable mirror to substantially conform 
to a calculated surface shape, 
wherein said first deformable mirror, said second deform- 
able mirror and said signal processing and control sys- 
35 tem together provide a large amplitude light modulation 
range to provide an actively shaped optical beam, 
wherein said large amplitude modulation range is charac- 
terized by 

40 


wherein 

45 

\dh\ 

I dr\ 
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is a focal gradient of at least one of said first deformable 
mirror and said second deformable mirror at a point of a 
respective reflecting surface thereof, 

wherein h is a surface deflection of a corresponding one of 
55 said first deformable mirror and said second deformable 
mirror at said point, and 

wherein 0 S is a solid angle subtended by a disk of radius 5 
and the distance z between said first deformable mirror 
and said second deformable mirror. 

60 2. .An active optical beam shaping system according to 

claim 1, wherein said large amplitude modulation range is 
greater than 1.2% of a frilly illuminated portion of said 
entrance beam of light. 

3. .An active optical beam shaping system according to 
65 claim 2, wherein said large amplitude modulation range is 
substantially an entire range to permit selectable amplitude 
modulation of portions of said entrance beam of light any- 
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where from substantially fully attenuated to substantially 
fully illuminated relative to fully illuminated portions of said 
entrance beam of light. 

4 . An active optical beam shaping system according to 
claim 1, further comprising a beam analyzer arranged to 
receive at least a portion of said exit beam of light reflected 
from said second deformable mirror, 

wherein said beam analyzer is configured to communicate 
with said signal processing and control system to pro- 
vide feedback control of at least one of said first and 
second deformable mirrors. 

5. An active optical beam shaping system according to 
claim 1, wherein said signal processing and control system is 
further configured to calculate said surface shape of at least 
one of said first and second deformable mirrors based on 
reverse ray tracing calculations. 

6. An active optical beam shaping system according to 
claim 1, wherein said signal processing and control system is 
further configured to calculate said surface shape of at least 
one of said first and second deformable mirrors based on a 
solution to Monge-Ampere equations corresponding to said 
active optical beam shaping system. 

7. An optical transmitter, comprising: 

an optical source; and 

an optical modulator arranged in an optical path of said 
optical source, 

wherein said optical modulator comprises: 

a first deformable mirror arranged to at least partially 
intercept an entrance beam of light from said optical 
source and to provide a first reflected beam of light; 
a second deformable mirror arranged to at least partially 
intercept said first reflected beam of light from said 
first deformable mirror and to provide a second 
reflected beam of light; and 
a signal processing and control system configured to 
communicate with said first and second deformable 
mirrors, 

wherein said signal processing and control system is con- 
figured to provide control signals to said first deformable 
mirror and said second deformable mirror to provide at 
least a portion of said second reflected beam of light as 
an output beam having at least one of a selectable ampli- 
tude profile or selectable phase profile, 

wherein said first deformable mirror, said second deform- 
able mirror and said signal processing and control sys- 
tem together provide a large amplitude light modulation 
range to provide an actively shaped optical beam, 

wherein said large amplitude modulation range is charac- 
terized by 


wherein 

I 

lari 

is a focal gradient of at least one of said first deformable 
mirror and said second deformable mirror at a point of a 
respective reflecting surface thereof, 

wherein h is a surface deflection of a corresponding one of 
said first deformable mirror and said second deformable 
mirror at said point, and 
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wherein 0^ is a solid angle subtended by a disk of radius 8 
and the distance z between said first deformable mirror 
and said second deformable mirror. 

8. An optical transmitter according to claim 7, wherein said 
5 signal processing and control system is configured to provide 

control signals to said first deformable mirror and said second 
deformable mirror such that said output beam lias both a 
selectable amplitude profile and selectable phase profile. 

9. An optical transmitter according to claim 8, wherein said 
to signal processing and control system is configured to provide 

control signals to said first deformable mirror and said second 
deformable mirror such that said output beam of light com- 
prises soliton pulses. 

10. An optical receiver, comprising: 

15 an optical detector; and 

an optical filter arranged in an optical path of said optical 
detector, 

wherein said optical filter comprises: 
a first deformable mirror arranged to at least partially inter- 
20 cept an entrance beam of light being detected to provide 

a first reflected beam of light; 
a second deformable mirror arranged to at least partially 
intercept said first reflected beam of light from said first 
deformable mirror and to provide a second reflected 
25 beam of light; and 

a signal processing and control system configured to com- 
municate with said first and second deformable mirrors, 
wherein said signal processing and control system is con- 
figured to provide control signals to said first deformable 
30 mirror and said second deformable mirror to provide at 
least a portion of said second reflected beam of light as 
an output beam directed to said optical detector having at 
least one of a selectable amplitude profile or selectable 
phase profile, 

35 wherein said first deformable mirror, said second deform- 
able mirror and said signal processing and control sys- 
tem together provide a large amplitude light modulation 
range to provide an actively shaped optical beam, 
wherein said large amplitude modulation range is charac- 
40 terized by 



45 

wherein 


\dh\ 


is a focal gradient of at least one of said first deformable 
mirror and said second deformable mirror at a point of a 
55 respective reflecting surface thereof, 

wherein h is a surface deflection of a corresponding one of 
said first deformable mirror and said second deformable 
mirror at said point, and 

wherein 0^ is a solid angle subtended by a disk of radius 8 
60 and the distance z between said first deformable mirror 
and said second deformable mirror. 

11. An optical communication system, comprising: 
an optical transmitter; 

an optical transmission path optically coupled to said opti- 
65 cal transmitter; and 

an optical receiver optically coupled to said optical trans- 
mission path. 
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wherein at least one of said optical transmitter, said optical 
transmission path, or said optical receiver comprises an 
active optical beam shaper, and 
wherein said active optical beam shaper comprises: 

a first deformable mirror arranged to at least partially 5 
intercept an entrance beam of light to provide a first 
reflected beam of light; 

a second deformable mirror arranged to at least partially 
intercept said first reflected beam of light from said 
first deformable mirror and to provide a second 
reflected beam of light; and 
a signal processing and control system configured to 
communicate with said first and second deformable 
mirrors, 

wherein said signal processing and control system is 
configured to provide control signals to said first 15 
deformable mirror and said second deformable mirror 
to provide at least a portion of said second reflected 
beam of light as an output beam having at least one of 
a selectable amplitude profile or selectable phase pro- 
file, 20 

wherein said first deformable mirror, said second 
deformable mirror and said signal processing and 
control system together provide a large amplitude 
light modulation range to provide an actively shaped 
optical beam, 25 

wherein said large amplitude modulation range is charac- 
terized by 


ah 

dr 


* 0 ( 9 ,) 


30 


wherein 


35 


ah 

dr 
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15. An optical telescope, comprising: 
a light collection and magnification system; 
an active optical beam shaper arranged in an optical path of 
light collected and magnified by said light collection and 
magnification system; and 

an optical detection system arranged to receive a beam of 
light from said active optical beam shaper, 
wherein said active optical beam shaper comprises: 
a first deformable mirror arranged to at least partially 
intercept an entrance beam of light to provide a first 
reflected beam of light; 

a second deformable mirror arranged to at least partially 
intercept said first reflected beam of light from said 
first deformable mirror and to provide a second 
reflected beam of light; and 
a signal processing and control system configured to 
communicate with said first and second deformable 
mirrors, 

wherein said signal processing and control system is con- 
figured to provide control signals to said first deformable 
mirror and said second deformable mirror to provide at 
least a portion of said second reflected beam of light as 
an output beam having at least one of a selectable ampli- 
tude profile or selectable phase profile, 
wherein said first deformable mirror, said second deform- 
able mirror and said signal processing and control sys- 
tem together provide a large amplitude light modulation 
range to provide an actively shaped optical beam, 
wherein said large amplitude modulation range is charac- 
terized by 


wherein 


is a focal gradient of at least one of said first deformable 
mirror and said second deformable mirror at a point of a 
respective reflecting surface thereof, 

wherein h is a surface deflection of a corresponding one of 
said first deformable mirror and said second deformable 
mirror at said point, and 

wherein 0^ is a solid angle subtended by a disk of radius 8 
and the distance z between said first deformable mirror 
and said second defonnable mirror. 

12. An optical communication system according to claim 
11, wherein said optical transmission path is at least one of an 
atmospheric or substantially a vacuum transmission path to 
provide a free-space optical transmission path, 

wherein said an optical transmitter is arranged to transmit 
a free-space optical signal, and 

wherein said optical receiver is arranged to receive a free- 
space optical signal such that said optical communica- 
tion system is a free-space optical communication sys- 
tem. 

13. An optical communication system according to claim 
11, wherein said optical transmission path comprises an opti- 
cal waveguide arranged to receive optical signals from said 
optical receiver and transmit said optical signals to said opti- 
cal receiver. 

14. An optical co mmu nication system according to claim 
13, wherein said optical waveguide is an optical fiber such 
that said optical communication system is a fiber-optic com- 
munication system. 


ia/!i 
I dr I 

is a focal gradient of at least one of said first deformable 
mirror and said second deformable mirror at a point of a 
respective reflecting surface thereof, 

wherein h is a surface deflection of a corresponding one of 
said first deformable mirror and said second deformable 
mirror at said point, and 

wherein 0 S is a solid angle subtended by a disk of radius 8 
and the distance z between said first deformable mirror 
and said second defonnable mirror. 

16. A solar furnace, comprising: 
a light collection system; 

an active optical beam shaper arranged in an optical path of 
light collected by said light collection system; and 
an optical focusing system arranged to receive a beam of 
light from said active optical beam shaper to be focused 
on an object to be heated, 

wherein said active optical beam shaper comprises: 
a first defonnable minor arranged to at least partially 
intercept an entrance beam of light to provide a first 
reflected beam of light; 

a second defonnable mirror arranged to at least partially 
intercept said first reflected beam of light from said 
first defonnable minor and to provide a second 
reflected beam of light; and 
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a signal processing and control system configured to 
communicate with said first and second deformable 
mirrors, 

wherein said signal processing and control system is con- 
figured to provide control signals to said first deformable 5 
mirror and said second deformable mirror to provide at 
least a portion of said second reflected beam of light as 
an output beam having at least one of a selectable ampli- 
tude profile or selectable phase profile, 
wherein said first deformable mirror, said second deform- 
able mirror and said signal processing and control sys- 
tem together provide a large amplitude light modulation 
range to provide an actively shaped optical beam, 

wherein said large amplitude modulation range is charac- 
terized by 


wherein 

I ^1 

I dr I 

is a focal gradient of at least one of said first deformable 
mirror and said second deformable mirror at a point of a 
respective reflecting surface thereof, 
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wherein h is a surface deflection of a corresponding one of 
said first deformable mirror and said second deformable 
mirror at said point, and 

wherein 0^ is a solid angle subtended by a disk of radius 8 
and the distance z between said first deformable mirror 
and said second deformable mirror. 

17. An optical pulse-shaping system, comprising: 
an optical pulse source; 

a first diffraction grating disposed in an optical path of said 
optical pulse source; 

an active optical beam shaper disposed in an optical path of 
light diffracted from said first diffraction grating; and 
a second diffraction grating disposed in an optical path of 
light output from said active optical beam shaper, 
wherein said active optical beam shaper comprises: 
a first deformable mirror arranged to at least partially 
intercept an entrance beam of light to provide a first 
reflected beam of light; 

a second deformable mirror arranged to at least partially 
intercept said first reflected beam of light from said 
20 first deformable mirror and to provide a second 

reflected beam of light; and 
a signal processing and control system configured to 
communicate with said first and second deformable 
mirrors, 

25 wherein said signal processing and control system is con- 
figured to provide control signals to said first deformable 
mirror and said second defonnable mirror to provide at 
least a portion of said second reflected beam of light as 
an output beam having at least one of a selectable arnpli- 
30 tude profile or selectable phase profile. 



